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ABSTRACT 


In  this  report  we  summarize  the  research  performed  under  the 
ONR  Grant  N00014-75-C-0293  during  the  calendar  years  1979  and  1980. 

We  describe  a  number  of  electromagnetic  scattering  problems  which 
we  have  solved  using  the  spectral  domain  methods.  Among  these  are 
the  frequency  selective  surfaces,  lealcy-wave  antennas  on  dielectric 
waveguides  and  arbitrarily  shaped  conducting  or  dielectric  scatterers 
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I .  INTRODUCTION 


This  research  effort  was  to  develop  new  approaches  to  solving  electro¬ 
magnetic  and  acoustic  scattering  problems  ^n  frequency  regimes  and  for 
geometrical  configurations,  for  which  the  conventional  numerical  or 
asymptotic  techniques  are  found  to  be  inadequate,  inefficient  and/or 
inaccurate.  Rather  than  employ  either  the  matrix  met. -ad;  or  the  ray 
techniques,  which  are  typically  valid  in  the  low-  and  h  gh-frequency 
regimes  respectively,  we  investigate  the  application  of  the  FFT  algorithm 
in  the  spectral  o-  transform  domain.  Two  different  variations  of  the 
method  have  been  studied.  The  first  of  these  employed  asymptotic 
solutions,  such  as  those  based  on  the  Physical  Optics  or  GTD  methods,  as 
initial  approximations  for  an  iterative  procedure  for  constructing  the 
solution  to  the  scattering  problem.  The  second  version  of  the  spectral 
technique  utilized  the  variational  principle  and  developed  a  procedure 
called  the  spectral-Galerkln  method.  The  end  result  of  the  application 
of  the  latter  method  is  a  matrix  equation  for  the  coefficients  of  the 
expansion  functions  used  to  repr  ,««nc  the  nknown  field.  Typically,  the 
size  cf  this  matrix  is  much  smaller  than  the  one  obtained  via  conventional 
procedures.  Our  objectives  were  to  investigate  these  spectral  domain 
techniques  in  great  detail  and  to  evaluate  the  scope  and  limitations  of 
the  two  approaches. 
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II.  SUMMARY  OF  TECHNICAL  ACCOMPLISHMENTS 


During  Che  last  grant  period,  we  have  carried  out  an  extensive 
investigation  of  the  spectral-Galerkin  and  spectral  iteration  techniques 
and  have  applied  them  to  the  problem  of  analyzing  a  number  of  electromagnetic 
radiation  and  scattering  problems.  We  have  applied  the  spectral  approach 
to  the  problem  of  scattering  by  frequency  selective  surfaces  (FSS)  which 
find  widespread  use  in  radomes,  reflector  antennas,  and  optical  filters. 

The  spectral  iteration  technique  has  been  found  useful  in  the  low  and  inter¬ 
mediate  frequency  ranges  where  the  cell  size  or  the  period  of  the  FSS,  which 
comprises  a  screen  with  periodic  perforations,  is  on  the  order  of  two 
wavelengths  or  less.  Beyond  this  limit,  the  spectral-iteration  technique 
is  more  efficient  as  it  avoids  matrix  inversion  altogether  and  derives 
the  solution  to  the  integral  equation  using  an  iterative  procedure.  Using 
these  two  methods,  we  have  successfully  analyzed  several  different  versions 
of  FSS  and  have  compared  the  results  with  theoretical  and  experimental 
data  published  elsewhere.  The  results  derived  with  the  spectral  approach 
have  not  only  been  found  to  be  accurate  and  efficient,  but  are  also  applicable 
in  a  frequency  range  which  is  considerably  wider  than  that  of  the  conventional 
method. 

The  spectral  iteration  approach  has  also  been  found  useful  for 
analyzing  other  stucturea,  such  as  conducting  and  dielectric  scatterers 
of  arbitrary  shape.  To-date,  only  a  preliminary  investigation  of  this 
problem  has  been  carried  out  but  the  results  appear  to  be  quite  encouraging. 

An  invited  paper  describing  the  spectral  technique  was  recently  presented 
at  the  Method  of  Moments  Workshop  held  in  St.  Cloud,  Florida,  under  the 
auspices  of  Rome  Air  Development  Center.  A  paper  describing  this 


presentation  is  appended  herewith  (Appendix  I) .  Other  manuscripts  accepted 
for  journal  publication  and  describing  the  work  carried  out  with  partial 
support  from  this  grant  are  also  attached. 
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APPENDIX  A 


SOLVING  ELECTROMAGNETIC  SCATTERING  PROBLEMS 
WITHOUT  MATRIX  INVERSION 

R.  Mittra,  C.  H.  Tsao  and  R.  Kastner 
Electromagnetics  Laboratory 
University  of  Illinois 
Urbana,  IL  61801 

Abstract  -  The  applications  of  the  Moment  Method  a  la  Harrington 
to  the  solution  of  electromagnetic  scattering  and  radiation  problems 
are  well  known,  and  the  method  has  revolutionized  the  way  boundary- 
value  problems  are  being  solved  today  on  modern  computers.  However, 
as  the  frequency  becomes  higher  and  the  body  size  becomes  comparable 
to  the  wavelength  of  the  incident  field,  the  CPU  time  on  the  computer 
becomes  large  and  the  storage  requirements  also  become  large  -  if  not 
prohibitive.  It  is  therefore  useful  to  look  for  alternative  approaches 
to  the  moment  method  for  attacking  the  radiation  and  scattering  prob¬ 
lems  in  the  so-called  resonance  region  and  above. 

In  this  paper,  we  introduce  an  iterative  technique  in  the  spectral 
domain  which  circumvents  the  limitations  of  the  moment  method  alluded 
to  above.  The  method  is  computationally  efficient  because  it  makes  ex¬ 
tensive  use  of  the  FFT  algorithm  to  perform  the  Fourier  transformation, 
which  is  an  Integral  part  of  the  spectral  domain  approach.  The  pro¬ 
cedure  also  has  the  uni  ,ue  feature  that  it  has  built-in  convergence  and 
accuracy  checks,  features  which  are  not  typically  found  in  other  methods. 
The  psper  Illustrates  the  application  of  the  spectral-iteration  tech¬ 
nique  using  scattering  from  periodic  structures  and  arbitrary  bodies  as 
examples. 


The  work  was  supported  by  the  Office  of  Naval  Research,  Contract 
N00014-75-C-0293. 


I .  INTRODUCTION 


The  purpose  of  this  paper  is  to  describe  an  approach  called  the 
spectral-iteration  technique  for  solving  electromagnetic  scattering 
problems  without  the  need  for  matrix  inversion.  The  method  is  espe¬ 
cially  suited  in  the  high  frequency  range  where  the  dimensions  of 
structure  are  large  compared  to  the  wavelength.  If  the  moment  method 
were  applied  in  this  range,  the  matrix  size  that  would  be  required  to 
handle  such  structures  would  be  large,  and  the  matrix  inversion  time 
as  well  as  the  storage  cost  would  be  prohibitive.  Also,  there  are 
geometries  such  as  grating  structures  which  we  will  be  discussing 
shortly,  for  which  no  asymptotic  solutions  are  available  because  the 
ray  solutions  based  on  GTD,  or  physical  optics  approximations,  are 
entirely  inadequate. 

Although  the  spectral-iteration  technique  has  recently  been 
applied  to  a  wide  class  of  problems,*  for  the  sake  of  illustrating 
the  principles  of  the  method  we  will  use  the  example  of  periodic 
structures  such  as  arrays  of  conducting  strips  or  periodically  per¬ 
forated  screens  which  can  be  either  free-standing  or  printed  on 
dielectric  substrates  (see  Figure  1).  These  gratings  have  frequency 
selective  properties,  ahd  find  many  applications  as  artificial  die¬ 
lectrics,  optical  and  quasi-optlcal  devices,  and  dichrolc  surfaces 
for  antenna  reflectors  and  radomes. 

Conventionally,  the  problem  of  electromagnetic  scattering  from 
these  periodic  structures  is  attacked  using  the  mode-matching  pro¬ 
cedure  employed  in  conjunction  with  the  method  of  moments.  A  de¬ 
scription  of  this  procedure  can  be  found  in  a  number  of  papers  on  the 
subject  by  Chen  [1],  Lee  [2],  and  McPhedran  and  Maystre  [3].  Though 
this  method  works  quite  well  in  the  low-frequency  region,  it  becomes 
prohibitively  costly  if  not  impractical  at  the  high  frequency  region 
where  the  aperture  size  is  one  to  two  wavelengths,  or  larger,  because 


For  a  bibliography  on  Spectral  Domain  Methods  refer  to  the 
Appendix. 
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Fleur*  1.  Freauenev  selective  surface. 


the  matrix  size  required  for  an  accurate  solution  becomes  prohibi¬ 
tively  large  and'  tt\£  numerical  computation  becomes  extremely  time- 
conr  -sing  and  costly.  As  mentioned  earlier,  the  high  frequency 
techniques,  e.g.t  GTD,  cannot  be  applied  to  circumve it  the  above 
difficulty  either,  because  the  complex  geometrical  configuration 
of  the  structure  does  not  lend  itself  to  the  ray  fo-malism  of  GTD. 

In  this  paper  we  introduce  a  new  technique  based  on  the  spectral 
domain  approach  which  provides  an  efficient  and  accurate  solution 
to  the  grating  problems  described  in  this  paper. 

As  a  first  step,  the  new  approach  begins  with  the  formulation 
of  the  problem  in  terms  of  an  integral  equation  in  the  transform 
domain.  The  standard  procedure  for  deriving  the  integral  equation 
for  the  unknown  aperture  field  (or  the  induced  current)  is  still 
followed;  however,  in  the  transform  domain  the  convolution  form  of 
the  Integral  equation  becomes  an  algebraic  one.  Furthermore,  be¬ 
cause  of  the  periodic  nature  of  the  structure,  the  transform  natur¬ 
ally  takes  the  form  of  DFT  (discrete  Fourier  transform)  which  can, 
in  turn,  be  efficiently  evaluated  using  the  FFT  (fast  Fourier  trans¬ 
form)  algorithm.  The  transformed  Integral  equation  is  subsequently 
solved,  using  an  Iterative  procedure,  simultaneously  for  the  aper¬ 
ture  field  and  the  Induced  current.  It  is  evident  that  the  method 
avoids  the  time-consuming  steps  of  evaluating  the  matrix  elements 
and  their  subsequent  inversion.  More  importantly,  the  problem  of 
storing  and  handling  over-sized  matrices  is  circumvented  even  at 
high  frequencies,  where  the  number  of  unknowns  can  exceed  the  figure 
2000.  An  added  feature  of  the  method  is  that  a  built-in  step  in  the 
iterative  procedure  provides  a  convenient  measure  for  the  boundary 
condition  check,  a  feature  not  readily  available  in  conventional 


problem.  The  above  authors  have  used  the  GTD  solution  as  the  zeroth 
order  approximate  solution  and  have  also  employed  an  Iterative  pro¬ 
cedure  to  generate  the  final  solution.  However,  to-date  this  pro¬ 
cedure  has  not  been  applied  to  the  grating  problems  being  considered 
in  this  paper. 

In  the  next  section,  we  present  the  formulation  of  the  periodic 
grating  problem.  In  section  III,  we  describe  the  iterative  procedure. 
In  section  IV  we  illustrate  the  application  of  the  technique  to  a 
number  of  practical  geometries.  Finally,  we  demonstrate  in  section  V 
that  the  approach  is  useful  for  a  class  of  closed-region  problems, 
e.g.,  waveguide  discontinuities.  A  brief  summary  of  the  paper  is 
Included  in  section  VI  and  some  conclusions  are  presented. 


II .  FORMULATION 

For  the  sake  of  illustrating  the  spectral  approach,  we  consider 
the  problem  of  a  uniform  plane  wave  scattered  from  a  free-standing 
periodically  perforated  conducting  screen  shown  in  Figure  2.  However, 
the  method  of  solution  is  easily  and  conveniently  extendable  to  the 
case  of  a  screen  on  a  dielectric  substrate. 


Due  to  the  periodicity  of  the  structure,  the  electric  field  on 
either  side  of  the  screen  can  be  expanded  in  terms  of  the  Floquet 
space  harmonics.  Using  the  e^(t  time  convention  (suppressed),  we  can 
write 
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for  z  <  0 
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where  E*  represents  the  incident  field,  (X+  ,  Y+  )  and  (X  ,  Y  ) 

am  ton  ran  ran 

are  the  reflection  and  the  transmission  coefficients  of  the  Floquet's 

harmonic  modes,  respectively  and  ^  ' s  stand  for  the  Floquet  mode 

mn 

functions,  given  by 

where 


v 
mn 

and 

V 

The  z-component  of  the  E-field  can  be  derived  from  (1)  using  the 
divergence  theorem. 

Enforcing  the  condition  that  the  tangential  field  is  continuous 
across  the  interface,  we  obtain 


^nn  "  exPU(umox  +  vmny)] 
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The  H-fleld  in  the  region  z  >  0  and  z  <  0  can  be  derived  from 
(1).  Evaluating  the  H-fleld  at  z  -  0-  and  z  -  0+,  subtracting  the 
expression  for  one  from  the  other,  and  making  use  of  the  requirement 
that  the  tangential  components  be  continuous  across  the  aperture, 
one  arrives  at  the  equation: 
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for  z  ■  0 
(2) 


I1  is  Che  incident  H- fie Id, 

J  is  Che  induced  current  on  Che  surface, 
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In  (2)  we  have  used  Che  nocacion  chat  for  a  function  f(r)  defined  on  Che 
z  ■  0  plane,  where  r  is  Che  position  vector  on  that  plane,  the  trunca¬ 
tion  operaCor  9  is  defined  by 


and 


0(f(r))  -  f(r)  for  r  on  Che  conducting  surface 

■  0  for  r  in  the  aperture 

§<f(r))  -  f(r)  -  (f (?)) 


The  obvious  identity  8(J)  -  J  and  that  z  *  [H(z  ■  0+)  -  H(z  ■  0  )]  ■  J 
have  also  been  used  in  deriving  (2). 


Unlike  the  integro-diff erential  equation  in  the  conventional  aethod, 
which  applies  only  in  the  aperture  (or  strip)  region,  (2)  is  valid  over 
the  entire  surface.  The  price  paid  for  extending  the  equation  to  the 
full  range  is  the  introduction  of  an  extra  unknown  J.  However,  as  we 
will  soon  see,  the  additional  unknown  J  can  be  solved  for  along  with  the 
aperture  field  using  the  iterative  procedure  discussed  in  the  next  section. 


III.  ITERATIVE  PROCEDURE 

The  suamatlon  Involved  in  (1)  and  (2)  can  be  readily  identified  as 
the  OPT  operation.  Let  F  be  the  operator  representing  the  DPT,  and  let 
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Ml 


represent  the  tangential  electric  field  in  the  transformed  domain. 

Identifying  (X  ,  Y  )  in  (2)  as  the  Fourier  coefficients  £  .  and 
mn  sn  t 

writing  G  for  the  matrix 


we  can  write  (2)  symbolically  as 

F(S  •  i£)  -  -H*  +  §(5)  (3) 

where  the  subscript  t  indicates  the  tangential  components,  and  it  is 
understood  that  all  the  quantities  are  evaluated  at  z  ■  0. 

If  the  Induced  current  were  available,  the  solution  for  could 
be  immediately  obtained  by  Invoking  (3)  and  by  using 

lt  -  8"1(F_1(-hJ  +  §(3)))  (4) 

In  practice,  however,  J  is  the  unknown  to  be  solved  for,  together 

with  Ef  and  hence  (4)  cannot  be  used  directly.  Instead  of  using  (4), 

s  recursive  relation  between  the  (n+l)th  approximate  solution  g(n+^ 

a(n)  c 

and  the  nth  approximation  E£  is  now  derived  and  the  two  unknowns 
Ec  and  J  are  solved  for  simultaneously  using  an  Iterative  procedure. 

To  derive  the  recursion  formula  for  E^n\  we  begin  with  (3), 
which  relates  and  E^n\  and  write 

§<J(n))  -  F(G  •  E^a))  +  H*  (5) 

Substituting  (S)  into  (4),  one  obtains 

?tn+1)  ■  +  8(F(5  •  l<n))  +  H*j))  16) 

Equation  (6)  is  the  desired  recursive  formula.  Before  inserting 
into  (3),  we  adjust  its  amplitude  by  multiplying  with  a  scale 
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factor  K,  computed  according  to  the  variation  expression 


K 


<E 


p(n) 


-Hi> 
_ t 


<i(D) 

<Et  » 


P(2 


EcCn))> 


(7) 


where 


<f ,  »>  ■/  f  •  i 

*  *  aperture  * 


da. 


Equation  (7)  is  obtained  by  applying  the  one-term  Calerkln's 
method  to  (3)  using  E^  as  the  testing  function.  It  is  apparent 
that  K  ■  1  when  E^  is  the  exact  solution.  K,  therefore,  also 
provides  an  indication  of  the  accuracy  of  the  nth  iterated  result 
in  a  weighted-average  sense. 


In  the  following  we  proceed  to  outline  an  iterative  procedure 
for  solving  (6): 

1.  Begin  with  an  initial  estimate  E^.  The  amplitude  of 
-(0)  c 

Et  is  to  be  properly  adjusted  using  the  scale  factor 
K  determined  from  (7). 

2.  Compute  1^,  the  discrete  Fourier  transform  of  E^. 

This  step  car  be  carried  out  efficiently  using  the 
FFT  algorithm* 

3.  Compute  G  •  E^  . 

4.  Obtain  the  OFT  of  G  *  1^  using  FFT. 

5.  Subtract  -Ht  from  the  result  obtained  from  step  4. 

This  gives  the  zeroth-order  approximate  solution  J^. 
Generally,  the  approximate  solution  for  J  obtained  in 
this  step  has  non-zero  values  extending  beyond  the  con¬ 
ducting  surface.  The  satisfaction  of  the  boundary  con¬ 
dition  for  the  induced  current  can  thus  be  verified  by 
checking  how  well  the  nth  approximation  for  the  current 
3^  la  confined  to  the  conducting  surface. 

6.  Add  -1*  to  3^  obtained  In  the  last  step,  and  taka  the 
Inverse  OFT  of  the  result  using  FFT. 
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7.  Multiply  g-1  by  Ch«  result  obtained  from  step  6, 
obtaining  E^A/ 

8.  Take  the  Inverse  transfora  of  1^  to  get  E^. 

The  exact  solution  for  should  have  zero  value  on 
the  conducting  surface.  This  criterion  serves  as  a 
boundary  condition  check  for  the  approximate  solu¬ 
tion  obtained  in  this  step. 

9.  Repeat  the  whole  procedure,  as  necessary,  using 

9(2^)  to  generate  the  next  higher-order  solu¬ 
tions  and  E^  until  convergence  is  achieved. 

In  the  following  two  sections  several  examples  are  presented  to 
illustrate  the  application  of  the  technique  described  above. 

IV.  SCATTERING  FROM  GRATINGS  AMD  GRIDS 

Let  us  consider  a  free-standing,  strip  grating  structure  illumi¬ 
nated  by  a  normally  Incident  uniform  plane  wave  as  shown  in  Figure  3. 
Let  the  incident  E-field  be  polarized  parallel  to  the  ed*es  of  the 
strip  (an  H-wave). 

The  formulation  for  this  problem  is  given  by  (2).  The  iterative 
procedure  discussed  in  section  Ill  is  applied  to  solve  for  the  tan¬ 
gential  aperture  E-fleld,  E(,  and  the  induced  current  density,  J. 

Figure  4a  shows  the  Incident  E-fleld  truncated  in  the  aperture, 
which  is  used  as  the  zeroth-order  approximation  for  Ec>  i.e., 

E ^  ■  0(1*).  The  derived  from  E^  is  shown  in  Figure  5a. 
has  significant  non-zero  values  extending  into  the  aperture  region. 

This  could  be  expected  because  of  the  crude  initial  estimate  made  for 
Figure  5b  shows  obtained  after  one  iteration.  Observe  the 
significant  improvement  achieved  with  Just  a  single  iteration  even 
though  the  zeroth-order  approximation  for  E^  was  rather  crude.  Higher- 
order  solutions  for  1^  and  obtained  via  further  iterations,  are 

shewn,  respectively,  la  Figures  4b  to  4d  and  Figures  5b  to  5d.  The 
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Magnitude  of  the  aperture  K-fleld  distribution  for  the  geometry 
shown  In  Figure  3,  obtained  via  four  consecutive  Iterations. 


Figure  5.  Magnitude  of  tlie  Induced  current  density  distribution  on  the 
strip  for  the  geometry  in  Figure  3,  obtained  via  four 
consecutive  Iterations. 


Magnitude  of  the  Induced  current  density  distribution  on 
strip  for  the  geometry  In  Figure  3,  obtained  via  four 
consecutive  iterations. 


rapid  convergence  and  the  accuracy,  which  is  verified  by  the  boundary 
condition  check  of  the  solutions,  are  well-demonstrated  in  these 
figures.  The  induced  current  density  J  also  shows  the  expected  edge 
behavior,  i.e.,  it  becomes  large  at  the  edges  as  it  should  for  the 
incident  H-wave. 

Figures  6  and  7  show  the  solutions  for  and  J  obtained  after 
four  iterations  when  the  gratings  are  illuminated  by  an  obliquely 
incident  plane  wave  with  an  incident  angle  0  ■  30°  and  with  the  H- 
fleld  polarized  parallel  to  the  edges  of  the  strip  (in  the  E-wave 
case).  Again,  the  boundary  conditions  are  satisfied  extremely  well 
by  the  results,  and  the  aperture  E-field  also  shows  the  expected  edge 
behavior  for  the  incident  E-wave. 

Next  we  consider  the  scattering  from  a  free-standing  conducting 
grid  illuminated  by  a  normally  Incident  plane  wave.  The  geometry  of 

the  problem  is  shown  in  Figure  8.  The  aperture  area  is  approximately 

2  2 
10A  whereas  the  cell  area  is  about  AAA  .  The  Initial  approximation 

for  E£  is  still  chosen  to  be  the  truncated  incident  field,  and  the 

dominant  component  of  the  tangential  aperture  E-field  is  shown  when 

the  incident  E-fleld  Is  polarized  in  the  y-direction. 

For  all  the  computations  in  this  section,  32  terms  in  the 
Floquet  expansion  functions  are  used  in  representing  the  unknown 
fields  along  each  of  the  two  dimensions.  This  leads  to  2*1  equiva¬ 
lent  unknowns  to  be  solved  for.  The  computation  time,  however,  re¬ 
quired  for  deriving  the  solution  is  quite  moderate  (5-6  secs,  of  CPU 
tlan  on  the  CDC  Cyber  175  System).  Clearly,  any  matrix  method  dealing 
with  such  a  large  number  of  unknowns  will  be  totally  impractical. 


=  0.7 


E-fleld  distribution  for  current  density  distribution 

a  free-standing  grating  on  the  strips  of  a  free- 

llluailnated  by  an  standinp  grating  illuminated 

obliquely  incident  E-  by  an  obliquely  incident  E- 


Figure  8.  Distribution  of  the  dominant  component  of  E-fieid  in  the  plane  of 
a  perforated  screen  with  rectangular  apertures.  (a)  E  -field 


I 


regions  of  dissimilar  dimensions.  An  open-region  type  example  is  a 
corrugated  surface  which  can  be  thought  of  as  a  junction  of  two 
regions,  viz.,  the  infinite  half-space  and  a  periodic  array  of  short- 
circuited  waveguides.  For  the  sake  of  si»aplicity,  a  closed-region 
type  problem  -  a  step  discontinuity  in  a  parallel-plate  waveguide  - 
is  considered  in  this  section.  The  geometry  is  shown  in  Figure  9. 

The  incident  field  is  a  TE  mode  wave.  The  formulation  of  this  prob¬ 
lem  can  be  found  in  the  literature.  The  integral  equation  is  given 
by 


b  b 

f  E  (x' )K~(x,x' )dx'  -  [  E  (x' )K+ (x,x* )dx*  -  2H1 

j  y  j  y  * 

o  o 

for  0<x<b,  z-0 


(8) 


til 


where 


and 


Ey  is  the  unknown  aperture  E-fleld, 

H*  is  the  incident  H-field, 
x 

K  (x,x')  ■  •  l  T  sin  ~  x  sin  ~  x' 

iw  Jl*  m  a  a 

00 

. .  1  r  _+  .  mir  mir  , 

K  (x,x  )  ■  — —  )  T  sin  —  x  sin  —  x 

jOJli _ y  _  m  n  k 


^J(k2  -  (f)V/2 


r  -< 

m 


((f)2  -  k¥/2 

v  a 


r  "J(k2  -  (y)V/2 


((f)2  -  k2)1/2 


if  k2  >  (f  )2 


if  k2  <  (— )2 

a 


if  k2  >  (f  )2 
0 


if  k2  <  (f)2 

D 


Note  that  (8)  is  defined  in  the  region  0  <  x  <  b.  To  apply  the 
iterative  technique,  ic  has  to  be  extended  to  the  full  range 
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0  <  x  <  a.  This  is  achieved  by  introducing  an  extra  unknown  func¬ 
tion  J(x),  and  th»  extended  equation  takes  the  fora 

a  a 


■ 

E  (x')K  (x,x')dx' 

y 

o 


|  a 

E  (x')K  (x,x’)dx'  +  8(-2H  )  +  §(J(x)) 

y  ^ 

o 


(9) 


for  0  <  x  <  a 


where  for  any  function  f(x) 


9<f(x))  -  f (x)  if  0  <  x  <  b 

-0  if  b  <  x  <  a 


and 

8<f(x))  -  f(x)  -  0(f (x)) 


/— xt  \ 

A  recursion  formula  relating  the  (n+l)th  order  solution  1  to 

the  nth  solution  E^  can  be  derived  via  a  procedure  similar  to  that 


developed  in  section  III. 


The  formula  is 

a 


(n+1) 


K  - 


E<n)K+  +  9(-2Hi)  +  §( 
y  ^ 


{  E<"V  -  ' 


E(n)K+-  0(-2H1)) 
y  x 


(10) 


Equation  (10)  is  now  solved,  using  an  iterative  procedure  similar 
to  that  developed  in  section  III.  The  Integrations  in  (10)  can  again 
be  carried  out  using  the  FFT  because  of  the  characteristic  nature  of 
the  kernels.  Figure  10  shows  the  aperture  E-field  distribution  at  the 
discontinuity.  The  initial  approximation  for  E^  is  taken  to  be  the 
Incident  field  truncated  in  the  aperture. 


The  result  is  obtained  in  three  iterations  with  32  expansion 
functions  used  in  representing  the  unknown  field.  We  note  that  the 
boundary  condition  on  E^,  viz.  ,  that  it  vanishes  at  the  edges,  is 
satisfied  by  the  Iterated  solution. 
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ield  distribution  at  the  step-discontinuity 


VI.  AN  EXTENSION  OF  THE  SPECTRAL- ITERATION  TECHNIQUE  FOR  ANALYZING 

SCATTERING  FROM  BODIES  OF  ARBITRARY  SHAPE 

In  the  past,  the  application  of  the  spectral  technique  has  been 
largely  restricted  to  scatterers  with  two-dimensional  geometries  and 
with  planar  facets. 

A  new  extension  of  the  method,  outlined  below,  allows  one  to 
analyze  the  problem  of  scattering  from  bodies  with  arbitrary  shape 
and  size.  The  method  has  wide  application  since  it  is  capable  of 
handling  body  shapes  which  are  not  conveniently  treated  either  with 
the  1  ow  or  the  high-frequency  techniques,  an  example  being  a  long 
and  narrow  cylinder  of  finite  length.  Moreover,  the  range  of  appli¬ 
cation  of  the  spectral-iteration  approach  encompasses  both  the  low- 
frequency  region,  where  it  can  provide  faster  and  more  efficient 
solution  than  the  conventional  moment  methods,  and  to  the  high- 
frequency  range  where  it  yields  both  the  surface  current  and  the 
near  fields.  In  addition,  the  method  has  a  built-in  boundary  con¬ 
dition  check,  a  unique  feature  not  present  in  asymptotic  techniques, 
e.g.,  the  ray  methods. 

The  first  step  in  the  spectral-iteration  method  is  to  repre¬ 
sent  the  scatterer  via  an  array  of  flat  slices.  The  current  on  each 
such  slice  is  sampled  by  a  thin  ring  on  the  left  side  of  the  slice 
as  shown  in  Figure  11.  The  iteration  begins  with  an  initial  guess 
for  the  current  distribution  on  the  scatterer  on  the  basis  of,  say, 
the  physical  optics  (PO)  approximation  or  the  GTD.  Next,  we  evalu¬ 
ate  the  field  at  the  first  plane  which  is  farthest  to  the  left  of  the 
scatterer.  A  process  of  improving  the  assumption  for  the  current  then 
starts.  It  is  done  one  slice  at  a  time  from  left  to  right,  as  explained 
below.  Once  the  whole  body  has  been  scanned  in  this  way,  the  first 
iteration  is  over,  and  the  process  is  repeated  until  a  convergence  is 
attained  and  the  boundary  conditions  are  satisfied. 

The  manner  in  which  the  current  rings  are  updated  is  described  as 
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follows.  Consider  Che  pth  slice,  one  which  is  bounded  by  che  pth 

plene  on  Che  left  and  Che  p+lth  plane  on  Che  right,  and  contains 

the  thin  current  ring  J  at  its  left  imedlately  to  the  right  of 

P 

Che  pth  plane  (Figure  11).  At  the  pth-plane  the  scattered  field 

••  ^  — 

consists  of  two  coaponents,  viz.,  and  E^.  is  the  aggregate 
of  the  contribution  from  all  of  the  sources  to  the  left  of  the  pth 
plane.  Itj  spectrum  is  thus  propagating  to  the  right  and,  conse¬ 
quently,  the  radiation  condition  dictates  the  choice  of  the  propa¬ 
gator  e  ^zz  for  any  of  its  plane-wave  spectral  components.  The 

— 

second  contributor  to  the  scattered  field  is  E  ,  which  is  the  con- 

P 

trlbutlon  of  all  the  currents  to  the  right  of  the  pth  plane  and  is 
thus  propagating  to  the  left.  We  assume  that  in  the  process  of 
scanning  the  body  from  the  left  to  right  all  the  currents  to  the 
left  of  the  pth  plane  have  been  updated,  implying  the  E~  has  been 
updated  as  well. 

We  next  proceed  to  update  the  current  in  the  pth  slice,  i.e., 

J  .  To  do  this,  we  look  at  the  plane  immediately  to  the  right  of 
Jp.  This  plane  is  shown  by  the  dashed  line  in  Figure  S.l.  can 

be  transformed  to  this  plane  simply  by  adding  the  contribution  of 
Jp.  In  the  spectral  domain  we  have  the  expression 

t*  -  E "  +  C  •  J  (11) 

P  P  P 

where  G  is  the  Green's  function  in  the  spectral  domain,  the  -  super¬ 
script  denotes  the  plane  immediately  to  the  left  of  J  ,  and  the  + 

P 

superscript  is  associated  with  the  plane  immediately  to  the  right 
of  J  .  Since  E  is  known,  we  are  now  able  to  use  the  assumed  J  in 

P  P  *+  P 

conjunction  with  (11)  to  compute  E  and  then  the  total  scattered 
field 

E  -  (E+  +  E+) 

P  P 

Next,  Inverse  transforming  E^  yields  the  scattered  electric  field  in 
the  spatial  domain  and  the  application  of  tha  boundary  condition 
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allows  one  to  replace  Che  tocal  scattered  field  E  inside  the  body 

inc  ^ 

by  -E  .A  Fourier  transform  is  then  taken  and  the  following  equa¬ 


tion  is  used  to  obtain  the  updated  E 


2  + 


E+  -  E+  (updated)  -  E  +  (12) 

P  P  P 

Finally,  the  updated  J  is  derived  from  (11)  and  the  updated  E+ . 

^  ■++  ^ 

Also,  E  , ,  is  obtained  from  E  via  the  equation 

p+i  p 

--  -+  -jv 

Vi"EP  6  (13) 

This  completes  the  operation  on  Che  pth  slice  and  we  move  to  the 

(p+l)th  slice  to  repeat  the  process  in  order  to  obtain  a  new  value  for 

J  . , •  We  continue  in  the  same  manner,  proceeding  to  the  successive 
p-rl 

slices  toward  the  right  until  we  are  finished  with  all  the  slices  and 
have  covered  the  entire  body.  The  end  result  of  this  series  of  steps 
is  a  complete,  updated  version  of  the  current  on  the  entire  scatterfer. 
Having  obtained  this,  the  first  iteration  is  completed,  and  the  whole 
process  can  be  repeated. 

The  iteration  process  is  continued  until  convergence  is  achieved, 
as  indicated  by  the  satisfaction  of  the  boundary  condition  on  and  in 
the  interior  of  the  scatterer. 


It  should  be  noted  that  since  the  two-dimensional  FFT  Is  used, 
even  for  a  three-dimensional  scatterer,  the  method  is  computationally 
efficient  and  its  storage  requirement  is  low.  Furthermore,  since  the 
boundary-condition  check  is  applied  at  each  stage  of  the  iteration, 
the  accuracy  of  the  final  (convergent)  result  is  guaranteed. 

Preliminary  studies  have  indicated  that  the  arbitrary  body  scheme, 
though  originally  conceived  in  connection  with  perfectly  conducting 
scatterers,  may  prove  useful  for  handling  dielectric  scatterers  as 
veil.  It  appears  that  the  method  may  be  generalizable  to  inhomoge¬ 
neous  dielectric  bodies. as  well,  although  more  work  remains  to  be  done 
to  determine  the  scope  and  limitations  of  this  approach. 


T  >' 
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A  SPECTRAL  DOMAIN  APPROACH  FOR  COMPUTING  THE  RADIATION 
CHARACTERISTICS  CF  A  LEAKY-WAVE  ANTENNA  FOR  MILLIMETER  WAVES 

R.  Mittra  and  R.  Kastner 

ELECTRICAL  ENGINEERING  DEPARTMENT 
UNIVERSITY  OF  ILLINOIS 
UriBANA,  ILLINOIS  61801 

ABSTRACT 

This  paper  deals  with  a  new  method  for  evaluating  the  complex 
propagation  constant  3  in  a  leaky-wave  structure  comprising  of  thin, 
metallic  rectangular  strips  etched  on  a  dielectric  rod  of  rectangular 
cross  section.  The  radiation  pattern  of  the  leaky  wave  antenna  can  be 
determined  once  3  is  known,  since  Re(B)  governs  the  direction  of  the  main 
beam  and  Im(3)  accounts  for  the  beamwidth  and  aperture  efficiency.  Tn 
addition,  the  knowledge  of  the  dependence  of  3  on  frequency  allows  one  to 
design  the  antenna  for  frequency-scanning  applications.  The  method 
employed  in  this  paper  is  based  on  the  spectral  domain  approach  which 
formulates  the  eigenvalue  problem  in  the  Fourier  transform  domain. 
Computed  results  are  shown  to  be  in  very  good  agreement  with  experimental 


measurements. 


1. 


I ntrod u ctlo n 


In  this  paper  v/e  describe  a  novel  method  for  eval  ating  the  complex 
propagation  constant  8  in  a  leaky-wave  structure  comprised  of  an  array  of 
thin,  metallic,  rectangular  strips  etched  on  a  dielectric  rod  of  rectangula' 
cross  section-  The  geometry  of  the  problem  is  shown  in  Figure  1.  This 
configuration  finds  useful  applications  as  a  frequency  scannable  antenna, 
particularly  at  millimeter  waves  where  the  antenna  can  be  conveniently 
integrated  with  dielectric-based  planar  integrated  circuits  [1,2,3]. 

The  complex  propagation  constant  8  along  a  leaky-wave  antenna 
determines  the  radiation  pattern  of  the  antenna.  Specifically,  Re(3)  governs 
the  direction  of  the  main  beam,  and  Im  (8)  accounts  for  the  beam  width  and 
aperture  efficiency.  In  addition,  when  the  dependence  of  6  on  frequency 
is  known,  it  is  possible  to  design  the  antenna  for  frequency-scanning 
applications. 

In  the  past,  the  determination  of  S  has  often  been  accomplished  via 
experimental  means,  mostly  by  near-field  probing  techniques.  Analytical 
evaluation  of  8  has  been  carried  out  for  structures  for  which  the  leaky 
wave  is  generated  from  guiding  structure  which  supports  a  fast  wave,  e.g., 
a  slotted  waveguide  (A).  For  such  a  structure,  8  differs  only  slightly 
from  the  guided-vave  propagation  constant  in  the  absence  of  the  slots  as 


sets  of  networks  are  required  for  every  hybrid  mode.  In  addition,  the  dis¬ 
continuities  must  be  incorporated  into  the  equivalent  network  by  lump-element 
representation,  which  by  themselves  may  require  the  solution  of  some  involved 
boundary-value  problems. 

In  this  paper  we  employ  a  method  based  on  the  spectral  approach  which 
formulates  the  eigenvalue  problem  in  the  Fourier  transform  domain.  A  des¬ 
cription  of  this  method  appears  below. 

2.  Formulation  of  the  Problem 

The  formulation  of  the  problem  is  based  on  the  spectral  domain 
approach  which  has  the  following  advantageous  feature.  The  Green's  function 
for  the  dielectric  substrate  region  is  conveniently  expressible  in  the  trans¬ 
form  or  spectral  domain  in  a  closed  form,  whereas  it  takes  a  complicated  form 
in  the  conventional  space  domain  approach. 

Referring  to  Figure  1,  let  us  consider  a  y-polarized  wave  traveling 

along  the  z-direction.  Because  of  the  periodic  nature  oc  he  geometry,  we 

ccn  express  the  fields  propagating  along  the  struc.ure  in  terms  of  Floquet 

space  harmonics  with  wave  numbers  8  (=8+  2mr/d),  where  8  is  the  complex 

n 

wave  number  we  are  seeking.  It  is  evident  that  only  a  finite  number  of  8n 
are  in  the  visible  region,  i.e.,  satisfy  the  criterion  |Re  Bnl  <  k,  where  k 
is  the  free-space  wave  number,  and  only  these  8n  contribute  to  the  leaky-wave 
radiation.  Typically,  these  leaky-wave  antennas  are  designed  such  chat 
n  ■  ±1.  As  alluded  to  earlier,  the  real  part  of  8+^  determines  the  direction 
of  the  main  beam  and  its  scanning  properties,  whereas  Che  imaginary  part 
determines  the  bearawidth  and  efficiency  of  the  antenna. 


36 


As  a  first  step  toward  attacking  this  problem,  we  replace  the 
dielectric  rod  with  a  slab  using  the  well-know  effective  dielectric 
constant  method  [6],  [7].  The  effective  dielectric  constant  is  given 


by 


'eff 


£ 


fk  , 

_x 

,k  j 


2 


where  k  is  the  free-space  wave  number  and  k^  is  the  fundamental  mode  with 
spatial  frequency  corresponding  to  the  cosine  variation  in  the  y-direction 
Inside  the  dielectric.  Next,  the  Green's  function  for  the  geometry  is 
constructed  in  the  spectral  domain,  because  the  expression  for  the  Green's 
function  in  the  transform  domain  is  considerably  simpler  than  the  corres¬ 
ponding  one  in  the  space  domain.  The  expression  for  G  ,  the  two-dimensional 
Fourier  transform  of  the  spatial-domain  Green's  function  for  the  E-field 
of  y-directed  point  source  located  at  a  distance  d  above  the  dielectric 
surface  (Figure  2),  is  formulated  in  a  fashion  similar  to  Collin's  [8]. 

It  is  given  by 


G(::,u, v)  * 

.  „  2.  <*,1  *  \2)e  -  (1s.l  -  W8  -J’S.l1*-'*1 

-j  V-  -  "  >  7Tk~r  ,  -i\,T 

<ba  +  K;1y  •  -  -  ^2>  -  '  <2) 


where 


k 


xl 


/l-u^-v*1 


x2 


£eff-u-v- 


37 


and  u  and  v  are  the  transform  variables  corresponding  to  the  y  and  2 
directions,  respectively.  These  are  normalized  with  respect  to  k;  hence, 
all  the  dimensions  are  expressed  in  terrcs  of  electric  length  (radians). 

In  our  case,  d»0.  The  tedious  Fourier  inversion  of  (2)  is  avoided 
by  keeping  the  rest  of  the  development  in  the  transform  domain.  The 
representation  for  the  electric  field  in  terms  of  the  current  distri¬ 
bution  on  a  strip  takes  then  a  simple  form  of  algebraic  product  (as 
opposed  to  convolution  in  the  space  domain) .  We  have 


E  ■  G(o,u,v)  •  J(o,u,v) 


(3) 


where  J  is  the  transform  of  the  assumed  current  distribution  on  the 
strip.  Sirfce  w  <  <  L,  it  is  reasonable  to  assume  a  cosine  variation 
in  the  y  direction  and  uniform  distribution  in  the  z  direction.  J 
is  then  proportional  to 


J 


(u,v) 


cos 


uL 

2 


<T>  -  (f>' 


(4) 


To  obtain  the  eigenvalue  equation,  ve  superimpose  the  electric 
fields  generated  by  the  periodic  array  of  scrips  and  apply  Che  boundary 
cc-.dition  chat  the  total  electric  field  on  a  representative  strip  is 
zero.  The  nC**  strip  is  characterized  by  an  amplitude  e  and  a 

displacement  of  tU  in  the  z-direction.  Therefore,  the  transform  of 
its  current  distribution  is 


n 


(5) 
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and  the  total  field  is 


E-C5[  eJ (v-B)nA  W 

n«— ** 

The  total  E-field  on  any  strip  is  now  equated  to  zero  in  the  Calerkin 
sense: 

/  /  E  •  J  dy  dz  -  0  (7) 

strips 

By  Parseval's  theorem,  (7)  is  also  expressible  by 


J  J  E  •  J  du  dv  -  0 


(9) 


Substituting  (6)  into  (8),  one  obtains 


r  t  ' .  .  1? .  .  r  j  (v-S )  n& 

j  /  C(u,v)  J2(u,v)  l  e 


du  dv  m  0 


(10) 


Utilizing  the  identity 


l  ej(/" 3 i  Hv-Hf] 


(U) 


ore  c!  the  integrals  is  eliminated  and  (9)  becomes 


l  I  [G(u,v)  32(u,v))v_5  +  2n^  du  -  0  (12) 

n— •  A 

Equation  (12  )  is  the  desired  eigenvalue  equation.  The  solution  of 
complex  S  satisfying  (12)  yields  the  desired  complex  wave  ounber  foe 
Che  leaky-wave  antenna. 
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3.  Results. 

Equation  (12)  requires  an  integration  along  the  real  axis  in  the 
u-plane  and  a  summation  over  ail  Floquet  modes.  Since  all  the  consti¬ 
tuents  of  (12)  are  symmetric  with  respect  to  real  u,  the  integration 
need  be  carried  out  only  along  half  of  the  u-axis.  Moreover,  there 
appears  to  be  little  or  no  contribution  to  the  expression  beyond  u  -  2 

and  [n|  -  2. 

The  solucton  for  S  is  obtained  via  a  search  procedure  which 

seeks  the  zeroes  of  Eq.  (12)  within  some  numerical  tolerance.  The  antenna 

under  consideration  had  29  strips  with  the  following  dimensions  at  the  center 

frequency  of  f  -  80  GHz  (see  Figure  1):  L  -  0.8X,  T  -  0.3F67X,  W  -  0.3387X, 

d  ■  0.6667X,  and  with  e  -  2.46.  The  computed  6  for  f  »  0.95fQ,  1.05fQ  and 

l.lf  was  -.31  -J . 03 ,  -.222  -j.05,  -.13  -J.03  and  -.05  -J.03  respectively, 
o 

These  solutions  account  for  main  beam  directions  of  108°,  103  ,  97.5  and  93 
ith  respect  to  the  z-axis,  compared  to  the  experimental  values  of  109°, 

103°  98°  ’.nd  92°.  These  values  of  3  correspond  to  a  backward  wave  (n  -  -1). 
The  real  part  of  the  fundamental  (n  -  o)  B  differs  very  little  from  that  of 
an  unloaded  dielectric  rod  and  can  be  predicted  by  approximate  methods  [1-3]. 
However,  the  imaginary  part  of  6  has  not  been  analytically  computed  before 
for  this  structure.  The  experimental  results,  obtained  by  near-field  probing 
techniques,  yield  Ia(3)  *  -.03  which  is  in  very  good  agreement  with  our 
theory.  As  seen  from  Figure  3,  the  measured  and  computed  beam  widths,  of 
about  5°,  also  agree  well  with  each  other.  The  behavior  of  the  experimental 
pattern  in  the  side-lobe  region  is  attributable  to  the  radiation  from  the 
pattern  feed  region  which  was  not  modeled  on  the  theoretical  calculations. 
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Gecaetry  for  the  Construction  of  the  Green's  Function. 


Wynn 
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A  SPECTRAL  DOMAIN  ANALYSIS 
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ABSTRACT 

A  septum  and  an  impedance  matching  poet  are  used  as  a  beam  equalizer 
In  an  open-ended  waveguide- feed  for  reflectors  used  in  satellite  communica¬ 
tions  systems.  The  performance  of  this  design  over  a  frequency  band  is 
evaluated  using  a  spectral  domain  approach.  The  computed  radiation 
patterns  in  the  E-  and  H-planas,  as  wall  as  the  results  for  the  impedance 
match,  are  presented  in  the  paper. 
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I.  Introduction 


Rectangular  waveguide  array  feeds  for  reflector  antennas  play  an 
important  role  in  the  design  of  satellite  communication  systems.  To  make 
the  radiation  pattern  more  symmetric  in  the  E-  and  H-planes  of  the  feed, 
a  beam  equalizer  is  needed.  The  design  used  in  this  case  is  a  septum  placed 
across  the  mouth  of  che  waveguide  such  chat  the  aperture  distribution  is 
reshaped  to  satisfy  the  new  boundary  conditions  imposed  by  the  septum. 
Consequently,  Che  H-plane  radiation  pattern  is  narrowed  to  approach  the 
E- plane  pattern,  thereby  achieving  the  beam  equalizing  effect.  However, 
the  introduction  of  such  a  septum  creates  an  Impedance  mismatch  problem  for 
che  feed.  To  alleviate  this  problem,  a  matching  post  is  placed  behind  the 
septum  so  that  che  reflection  back  into  the  waveguide  is  minimized. 

The  performance  of  this  design  over  the  desired  frequency  band  is 
evaluated  using  a  spectral  domain  approach,  or  more  specifically,  Galerkin's 
method  applied  in  the  spectral  domain  [1].  The  scattered  fields  on  both 
sides  of  the  beam  equalizer  are  represented  In  terms  of  their  Fourier 
transforms  or  spectra  which  can  be  related  to  the  Induced  surface  currents 
on  the  septum  and  the  post.  These  unknown  induced  currents  are  expanded  in 
terms  of  known  basis  functions  and  unknown  coefficients.  A  matrix  equation 
for  the  unknown  coefficients  is  derived  by  applying  the  boundary  conditions, 
and  che  moment  method  is  Chen  employed  in  the  spectral  domain  to  solve  for 
these  unknown  coefficients,  which  In  turn  give  che  answer  to  che  unknown 
scattered  fields.  The  scattered  fields  for  all  modes  obtained  in  chls 
manner  are  Chen  used  to  compute  che  reflection  and  transmission  coefficients 
for  each  node,  propagating  or  attenuated.  A  cacit  assumption  made  is  chac 
the  scattered  field  on  che  open-ended  side  of  che  waveguide  is  che  same  as 
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chat  in  an  infinitely  long  waveguide  containing  the  beam  equalizer.  In 
ocher  words,  the  truncation  effects  of  the  waveguide  are  ignored  in 
this  analysis.  The  transmission  coefficients  are  used  to  weight  Che 
radiation  field  due  to  each  mode  of  waveguide  and  Che  superimposed 
radiation  pattern  is  computed.  The  reflection  coefficients  are  used 
to  assess  the  impedance  matching  performance.  Numerical  results  indicate 
that  the  E-  and  H-plane  principally  polarized  patterns  are  equalized 
extremely  well  over  the  entire  frequency  band  of  operation  and  that  the 
impedance  matching  is  also  quite  satisfactory. 

11.  Analysis 

The  geometry  of  the  waveguide  with  septum  and  post  is  shown  in 
Figure  1.  Since  the  cross-section  of  the  post  is  very  small,  the  post 
is  modeled  as  a  narrow  scrip  to  simplify  the  analysis.  The  incident 
field  is  propagating  in  the  z-direction  towards  the  post  as  shown 
schematically  in  Figure  1.  There  are  surface  currents  induced  on  the 
septum  and  Che  post  due  to  the  incident  field.  The  scattered  fields 
radiated  by  these  induced  surface  currents  then  propagate  in  both  the 
z-direction  and  the  -z-direction,  giving  rise  to  the  transmitted  and 
the  reflected  waves,  respectively.  In  the  following  analysis,  Che 
truncation  effects  of  the  waveguide  at  zm 0  are  ignored,  as  though  the 
post  and  septum  were  located  in  an  infinite  guide. 

The  incident  field  in  the  waveguide  can  be  expressed  in  terms  of 
TE  and  TM  modes  in  the  usual  manner: 

TE  modes: 

nm 
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u2  ue 


z~_  •  uu/a  ; 

1  £.  HQ 


2-pu  -  8  /(ue) 

CJi  ran 


The  scaecered  fields  can  be  expressed  in  terms  of  their  Fourier 
spectra,  which  are  In  turn  related  to  the  Fourier  spectra  of  the  induced 
surface  currents.  The  Fourier  spectra  of  the  induced  surface  currents 
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are  Chen  solved  for  by  Che  moment  mechod  applied  in  Che  transformed, 
i.e.,  the  spectral  domain.  Specifically,  Galerkin's  method  is  used 
in  the  present  analysis  —  che  same  basis  functions  are  used  as  testing 
functions  in  the  moment  method.  Upon  solving  the  spectra  of  the 
Induced  surface  currents,  Che  scattered  fields  can  be  obtained  in 
a  straightforward  manner.  The  analytic  details  follow. 

Cc  rresponding  co  each  incident  mode,  the  scattered  fields  E3  and 

x 

3 

can  be  represented  in  the  following  form: 


E3  (b/2  <  y  <  b) 
E3  (0  <  y  <  b/2) 

c.X 


•  cos  (—  x)  /  f  (a) 

A  H 


sinfy  (b-y)  ]| 
n 


sin  y  y 


exp(-jaz)da. 


(2) 


[e3  (b/2  <  y  <  b; 


E*  (0  <  y  <  b/2) 

*»•  — « 

where 


sin 


/  «*<*> 


sin[y  (b-v)] 
n 


s  in  y  y 
n 


exp(-jaz)doi 


1/2 


Y  -  [k2-  (^-)Z  -  a2] 

Q  A 


O) 


and  f^n),  gQ(:0  are  che  unknown  Fourier  spectra  to  be  determined.  The 
expression  for  E3  is  then  obtained  from  the  Maxwell's  equation  ?  *  E  *  0, 
giving 


E3  (b/2  <  y  <  b)*1 

„  f"cos(yn(b-y) ) 

•  sin (^7  x)  j  h  (a)  j 

4  *  Cl 

— oo 

t*  (0 <r< h/2) 

-cos  y_  y 

n 

49  s 

exp(-jaz)da 

(4) 
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wich 


h  (a) 
n 


(~)  f  (a)  +  ja  g  (a) 

an  n 


Mow  chat  we  have  E,  we  can  find  H  by  using  the  curl-of-E  Maxwell's  equation. 
However,  in  anticipation  of  relating  the  H~field  with  the  induced  surface 
currents  on  the  septum  and  the  strip,  only  the  x-  and  the  z-  components  of 
the  H-field  are  computed,  giving 


HS  (b/2  <  y  <  b) 
x 

|h*  (0<y<  b/2) 


T"~""  sin (—  x)  j  exp(-jaz) 
jam  a  1 


•  n_v„  +  J3  cos  t Y  (b**y)  J] 


[  [v  g  (a)  -  ja  h  (a)]  cos  y  y 
n  n  n  n 


H’  (b/2<  y  <  b) 


H  (0  <  y  <  b/2)| 

2  J 


x)  /  exp(-jaz) 


Lv  f  (a)  +  (~)  h  (■*)]  cos  ('r(b-y)] 

an  an  a 


( *y  f  (a)  -  (— )  h  (*)]  cos  y  y 

a  a  a  a  a 


(6) 


The  induced  currents  on  the  septum  and  the  strip  can  also  be  expressed 


in  terms  of  Fourier  spectra.  For  each  incident  mode,  we  have 


Jz  (x,y  -  b/2,  z) 


Jx  (x,y  -  b/2,  z) 


r*  ™ 

sin(—  x)  /  j  (a)  exp(-jaz)  da 


cos  (■—■  x)  /  j^Ca)  exp(-jaz)  da 


(7) 


where  jz(a)  and  j^ta)  are  Fourier  spectra  of  the  induced  surface  currents 
to  be  related  to  the  spectra  of  the  scattered  H-field.  This  relationship 
is  obtained  by  enforcing  boundary  conditions  on  the  septum  and  the  strip 
and  can  be  written  as 


(b/2  *  =  )>  ♦ 


<!  ■  £  <‘/2 


",  (b/2 


+  e))  -  Hz  (y  -  ^  (b/2 


£)) 


-  e)) 


(8) 


where  e  is  a  positive  quantity. 

Substitution  of  (6)  and  (7)  into  (8)  leads  to  the  following  algebraic 
equations : 


C 


_  y  b 

!  7“—  cos  (-T-)  [y_  g  (a)  -  jah  (a)] 
jwu  2  an  a 


Y  b 

~-cos  (-2-)  [  Y  f  (a)  *  (^)  h  (a)] 
jiuu  2  n  n  an 


jz(a) 


!  j  (a) 
I  X 


(9) 
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Substlcucing  bo(...  fro.  (5,  into  (9)  and  -nlpwlacldg  th,  r.suitl„g 
leads  Co  th.  following  nacrlx  «,o.clon  for  f^c)  and  ^(c  ln  cerms  of  £[], 
transform  domain  currants  j  (a)  and  j  (a) . 


-Ja  (S)  k2  -  3LZ* j 

a  'y 


^  -  a2  ja  (— ) 


r'nH 


8n(a>j 


^UYn 
Y  b 

2cos  (-T-) 


Va) 


JJa) 


An  Inversion  of  (IQ)  giv«3 


f„(o) 


-ia(HI)  k2.nV  .  : 


2  Y  b  ! 

2k  Y  cos(-f-)  ; 

*  2  '  t2  -  a2  j,<*S> 


V*! 

V.  w* 


Fourier  transforms  of  the  currents  J#(«>  and  are  known,  we 

can  obtain  f^o)  and  grfa)  from  (u).  Subsequent!,.  (2),  (4)  and  (6)  can 

derive  the  scattered  fields  by  substituting  for  f^a)  and  g  (o) 

in  those  expresaions.  The  raining  task  is  to  solve  for  th.  transform" 

currents  j,(«)  and  J^o).  using  the  Galerkin 's  method  applied  in  the 
transform  domain. 

First  we  express  the  corresponding  snace 

8  space  domain  currents  in  terms  of 

a  linear  combination  of  a  set  of  .uir,Kia,  v  a  , 

suitable  basis  functions  with  unknown 

coefficients. 
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0 


qj (z)  exp  (jaz)  dz,  i  -  1,  2 . 


"  2t 


-c 


(14) 


Replacing  (13)  in  (11),  one  obtains  the  expressions  for  the  transform 

domain  functions  f  (a)  and  G  (a) : 

n  n 


fn(a) 


gn(a) 


~1<»u 


2k2Y 


cos 


(Ip 


(k2  -  [  \  A  p  (a)|  -  f  \  B  Q  (a)' 

U-0  1  1  j  a  [l-l  4  4 


j  a  [  1  AiPi(a)J  +  (k2‘a2) 


r  L 

I  B  Q  (a) 

2.-1 


J 

(15) 


Substituting  these  expressions  in  (2) ,  (4)  end  (6)  we  obtain  the  final 
results  for  the  scattered  fields  in  terms  of  the  unknown  coefficients 
A’s  and  3's.  In  order  to  find  these  unknown  coefficients,  we  enforce  the 
boundary  condition  for  the  electric  field  on  the  septum  and  the  strip, 
which  requires  the  total  tangential  component  of  the  electric  field  to  be 
zero,  i.e., 


•E 


le‘J 


at  y  •- 


(16) 


Using  (1),  (2),  and  (15)  in  (16)  we  derive  the  following  equations  for  A’s 
and  B's. 

[k2  -  (~)^1  l  A.  f  ~  tan  j  rnb  P  (o)  exp(-jaz)  da 

L  *  J  i-0  —  y  n  ^  2  1 

L  r*  y  b 

•lr  I  B  —  tan  (-p)  Q  (a)  exp(-jaz)  da 

*  Z«1  1  Ya  L  L 
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^  -2  — 
r  1  a  8 


- 22  JL_  «  sin (—■)  exp  (+  jS  z),  for  TM  modes. 

'jju  ^4  run  L  tin 

c 


2~  ~  h  sin(~)  exp  (+  J3  z),  for  TE  modes, 

b  ,  2  na  2  nm 


Y  b 


j  —  r  *  /  —  tan  (-2-)  P.  (a)  exp  (-jaz)  da 

*  1-0  1  -  Tn  2  1 


L  «  ..2  2.  y  b 

+  l  B  /  — =2-A  Cjm  (.__)Q^(a)  exp  (-Jaz)  da 

l»l  -<*  Yn 


r  2 

+  ^ e  sin  (2p)  exp  (+  J3  z),  for  TM  modes 


w  y  am 


nm 


V. 


,  for  TE  modes 


Mow  multiply  both  sides  of  (17)  by  the  basis  functions  p i(z)  and  qj 
integrate  over  z  to  obtain  the  following  equations. 

1  -  ,  Y  b  * 

(k2  -2-5-)  r  A  /  ~  tan  (-2-)  P  (a)  P  ,  (a)  da 

a  i-0  — »  Yn 

L  -  Y  b 

■  ]  7  I  S.  /  T"  tan  (-§-)  Q  (a)  ?*,(»)  dc 
*  i-1  *  —  Yn  *  1  1 


(17) 


(z)  and 
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f-»< T>3 

i  UU 


*— SS  ?  .  .in? 

>u  ,  2  nm  2 

k 

c 


k  ,  .air 

2(T>  ~2  hm.*l0T 

kc 


P*i'  (3n»>  j 

Pf  «W  ! 


rv  <»«> 

V  <«,»> 


,  i'-O, for  TM  modes 


,  i'-O, for  TE  modes 


The  upper  quaneleies  vichin  Che  square  bracket  in  (18)  are  associated 
with  positive  z  Incident  nodes,  and  the  lower  quantities  are  associated 
with  negative  •  incident  modes.  Both  cases  are  included  because  we  are 
lncensced  in  the  transmission  coefficients  as  well  as  the  reflection 
coefficients.  The  integrals  are  evaluated  by  numerical  integration,  and 
the  resulting  system  of  linear  equations  is  solved  as  usual  for  the 
unknowns  A ' s  and  B's  by  matrix  .aversion. 
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The  choice  of  basis  functions,  i.e.,  che  {/sand  q's,  is  based  on 


previous  experience  and  che  consideration  of  the  behavior  of  the 
currants  at  che  edges  of  che  septum  and  the  scrip.  Only  two  terms  for 
che  expansion  of  che  currents  on  the  septum  are  retained,  which  are 
believed  to  be  adequate  for  this  analysis.  The  basis  functions  and 
their  Fourier  transforms  are  shown  in  Table  1. 

Let  us  now  consider  che  case  of  interest,  namely,  one  in  which  che 
incident  fields  are  TEq^  modes.  We  have  n*0,  and  che  following  relations: 


Y 


2 

0 


2 

a 


(19) 

(20) 


When  n«0,  (18)  can  be  simplified  to  the  following 

Y0b 

2  -  tan 

Z  Ai  /  —  ~h  Ma)  P?,(a)  da 

i-0  1  —  Vj.  1  1 

2 


4  .  .mir. 

«  3in(T) 


0m 


<  V  | 


i’-0,  1,  2 


2  -  y  b 

S  Bi  /  Yo  C4n  (2)  Ql(0t)  Q*'(a)  da  "  °* 


(21) 


The  equations  for  A's  and  B's  are  uncoupled.  Therefore,  A's  and  B's  can 
be  solved  for  separately.  The  system  of  equations  for  B's  is  homogeneous. 
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which  leads  co  Che  conclusion  Chat  all  B's  are  zero.  Consequently ,  there 
are  no  components  of  the  current  in  the  z  direction  either  In  the  septum 
or  the  scrip.  In  the  following  development,  only  Che  system  of  equations 
for  A's  is  investigated.  Observe  that  the  solution  of  A's  involves  an 
inversion  of  a  3  x  3  matrix,  which  is  an  easy  task  for  the  computer.  How¬ 
ever,  the  evaluation  of  the  matrix  elements  Involves  nine  complex  integrals 
to  be  numerically  integrated.  Since  there  are  singularities  in  these  integrals, 
we  must  examine  the  integrands  carefully  to  make  sure  that  Che  numerical 
integration  Is  applied  correctly  to  give  accurate  results  in  spite  of  the 
singularities.  Therefore,  it  is  useful  to  write  the  expressions  appearing 
in  the  integrands  of  (21)  in  an  explicit  manner  as  follows: 

VS  •  4  '-2  <t>2i 

VS  Jo  (f >  [J’(t  *  d  *  t>  1 

Vi  Jo  <¥>  t-<¥)2>  *•■“>  ho<|  +  i  *  f>i 

•  J  J!  (f'*  exp  t-(T)2]  exp  Cjo(2  +  d  +  f }  1 

2 

V$  -  -J  exp  C-Ja(2  +  d  +  2)] 


Hoce  chat  P^T?|,  ^2?1’  eCc'  are  noc  interesC  here  because  Che  corresponding 

integrands  are  odd  and  therefore  the  integrals  involving  them  are  identically 
zero.  Now  observe  chat  the  integrand  in  (21)  contains  a  simple  pole  located 
at  a  -  6  on  Che  path  of  integration  for  0  <  a  <  »  ,  provided  that  the 
wave  number  satisfies  the  following  condition: 


ir  ,  .  it 

b <k<  3b 


This  is  che  case  when  the  incident  field  is  Che  dominant,  propagating  mode 
in  Che  waveguide.  Introduction  of  some  loss  in  the  medium  clarifies  the 
position  of  the  poles  along  the  integration  path,  as  shown  in  Fig.  2.  Each 


of  che  integrals  in  (21)  can  be  written  in  Che  following  form 


f  F(a)da  -  [  8(51 


F(a)da  + 


F(a)da  +  [  F(a)da  +  w  j  Res  (-6m) 

•3  J  3  Ui 

01  301 


r23oi 

Tj  Res  (Sq^)  *  j  F(a)da  + 

J  —a* 


F(a)da 


F(a)da  + 


F(a)da  +  irj  Res  (~3ni)  -  irj  Res  (3rt1) 


F(a)do  + 


[F(-2B01  -a)  +  F(a)]  da 


+  j  *  [F(a)  +  F(23q1  -  a)]  da 


+  j  F(a)da  + 


Tj  Res  (-fl01)  -  Tj  Re,  (601) 
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In  the  above  equation,  a  bar  across  Che  integral  sign  means  principal  value 
integration.  The  method  of  foldover  as  shown  makes  the  new  integrand  remain 
bounded  at  Che  singularity  of  the  old  integrand;  hence,  this  new  integral 
is  easily  evaluated  by  the  numerical  method.  Also,  since  the  integrand 


goes  to  zero  rapidly  as  a  becomes  large,  the  integration  limit  »  can  be 
replaced  by  a  large  number,  e.g.,  14  3^.  It  should  be  noted  that  if  the 
foldover  mechod  was  not  used,  the  integral  could  still  be  evaluated  numer¬ 
ically,  in  some  cases,  but  the  integration  limit  08  must  be  replaced  by  a  much 
larger  number  because  of  the  heavy  tail  of  the  integrand,  and  doing  so  requires 
increased  computer  time  and  the  results  are  less  accurate.  The  residues  given 
in  (24)  are  easily  calculated  and  are  given  by  the  general  form 


Res(t801) 


b‘(s301> 


Pt  («)  Pr 


i,  i'  -  0,  1,  2 


with  Pj.  (a)  P^„  (a)  given  in  (22). 


(25) 


Having  discussed  the  evaluation  of  the  matrix  elements  in  detail,  we  can 

proceed  to  solve  the  matrix  equation  (21)  to  obtain  the  A's.  Having  found 

these  A's,  we  compute  fg(a)  from  (15)  and  chan  calculate  the  scattered  fields 

from  (2).  For  n*0,  both  E*  and  E*  •  0,  and  there  is  only  ES ,  which  can  be 

z  y  x 

written  as 


*1  <T  <  r  <  M 


E*  (0  *  7  «  7> 


P  -j«uo 


I  J  U  E  A  P  (a) 

2y.  cos  | ynb I  i-0 


sin  [yq  (b-y)l 


sin  Yq  y 


exp  (-jaz)  da 


(26) 
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where  yQ  - 


k>  ja  | 


f  j/a^-kz  if  k<  ja | 

V. 


Lee  us  evaluate  (26)  ac  y  *  —  , 


_s ,  b. 
Ex(y1) 


-jkbZ. 


2 

r  a  / 

i-0  1  - 


«  can 


YOb 


Y0b 


P^(a)  exp(-jaz)  da 


(27) 


Since  Che  scaceered  fields  on  both  sides  of  the  beam  equalizer  are  given  by 
(27) ,  we  can  compuce  Che  transmission  and  reflection  coefficients  by 
normalizing  these  scattered  fields  to  the  incident  field  given  in  (1). 

The  reflection  coefficient  R  can  be  expressed  as 


«p(jS0iZ) 


4  h 


01 


2 

l  A  [2*j  Z  Res .(c)  I 
i-0  1  ■-1,3,5,..  • 


(28) 


where  the  poles  are  given  by 


*  -801  "  *  ^  “  W  =  -;m>3  "  j  _k2 

and  the  residues  at  these  poles  are  given  by 


Resi(?1) 


4  P^-6^)  exp  (jzfl01) 

b2  <-801) 
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tn  >  3 


4  p, (jx_)  exp  (x  z) 
i  in _ m 


b  ( j  X  ) 

m 


with  x  ■  -k^~ 

m  b 

tx  ^ 

P0(jxm)  *^=  exp  [xm(c+d+|)]  exp  1 

W  *  f  «P  *0  (t  *■> 

P2(jV  ■  f  exp  <?T>  J1!  (f  xrn) 


and  IQ,  1^  are  modified  Bessel  functions  of  the  first  kind.  In  deriving  (23), 
the  integral  in  (27)  has  been  evaluated  by  the  residue  theorem  with  the 
contour  closed  in  the  upper  half  of  the  a-plane.  By  the  same  token,  the 
transmission  coefficient  T  can  also  be  obtained,  except  this  time  the  contour 
is  closed  in  the  lower  half  of  the  a-plane. 


-*exp(j301z) 


4h 


01 


2 

Z  A  [-2rrj  Z  Res  ('  )  J 

i-0  m-1,3,5,..,  1  m 


(29) 


where  the  poles  are  given 

r  a  3  a  XZTfi 

’i  oi  V 


by 


’a>3  " 
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and  che  residues  at  these  poles  are  given  by 


Res^;^)  - 


4  Pi(B01)  exp  C-jS012) 


b2S 


01 


4  p,MxJ  exp  (-x  2) 
Re3i(Cm)  "  Res<(“:)xJ  -  - - - ^ - - — 


i'  n 


m  >  3 


b‘WxJ 

in 


with  x 


m 


/S 


-k* 


CX. 


P0  f'jxm)  *  exp  t-Xm(c+d+  exP  f(~r^ 


P1  <*V  *  t  exp  (-'T')  :0  (f  xm} 

P2  ('jxm)  *  ¥  axp  (~xm  f}  ["jIl  (f:<mn 

and  IQ,  are  modified  Bessel  funccions  of  che  first  kind.  Numerical 
results  Indicate  chat  minimum  reflecti~r  and  maximum  transmission  can  be 
achieved  if  che  separation  between  the  septum  and  che  scrip  is  0.1 
wavelength  for  the  gi«en  dimensions  in  Fig.  1. 

The  radiation  pattern  of  the  waveguide  fe^d  with  beam  equalizer 
can  now  be  computed  in  the  following  manner.  First,  the  aperture  field 
distribution  in  the  plane  containing  the  waveguide  mouth  is  estimated  by 
a  superposition  of  che  waveguide  mode  field  at  this  plane  with  each  mode 
being  weighted  by  the  corresponding  transmission  coefficients  for  che  modes 
Since  the  transmission  coefficients  have  been  evaluated  at  a  different 
reference  plane,  it  is  necessary  to  refer  these  transmission  coefficient; 
back  to  che  aperture  plane.  This  is  done  by  multiplying  the  transmission 
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coefficients  by  an  appropriate  correcting  factor,  which  is  a  phase  factor 
for  a  propagating  mode  and  is  an  exponential  factor  for  an  attenuated 
mode.  The  radiation  pattern  is  then  obtained  by  the  familiar  Fourier 
transform  relation  between  the  far  field  and  the  aperture  field.  For  TE 

m 

modes,  the  far  fields  are  given  by  the  following  expressions  [2]: 

1/2  ,  .  n2  .  ,  3  3 

_  /j,  (fab)  sine  r,  ,  nm  „  .  „  ,,  nm  ,,, 

e3  ■  -  (-)  — T  z -  ^  +  "IT" cos9  +  R  ^  cos9)j 

2X  r  k 

nm 


[  (.—  sino)2  -  (if  coso)2]  *  Ce ,d>) 
a  o  am 


1/2  2 

^  m  _  “  (^ab)  sin6  sinO  coso 

*  £  2  X3  r 


3  3 

•  [cos3  +  —r—  +  R  (COS0 - -^)  ]  "•>  (9  , o) 

<c  nm 


(30) 


with 


•nm(9*5) 

nm 


- 

sin(-^  sin3  coso  +  ~) 

A  Z 

sin (-~  sin6  sino  +  nj-) 

2  2 

sine  coso)  -  (r~) 

.  2  2 
(—  sine  sino)  -  (— ) 

_  w 

.  exp  {  -j  [kr  -  —  sin9  (a  cosp  +  b  sino)  -  (n+nrt-1)  —  ] 


} 


where  (r,9,0)  are  the  conventional  right-handed  spherical  coordinates.  The 
total  radiation  pattern  is  then  obtained  by  a  superposition  of  these 
Individual  mode  patterns  with  the  appropriate  transmission  coefficients 
referenced  at  the  aperture  plane.  Computed  patterns  and  numerical  results 


are  presented  in  the  next  section. 


Ill •  Computed  Results 

The  performance  of  the  square  waveguide  feed  with  beam  equalizer 
whose  dimensions  are  shown  in  Fig.  1  with  d  **  0.1X  is  evaluated  over  a 
frequency  band  of  operation.  The  value  for  c  has  been  chosen  experi¬ 
mentally  to  be  0.25X.  The  computed  radiation  patterns  in  the  E-  and 
H-planes  are  presented  in  Figure  3.  The  beam  equalization  is  quite 
satisfactory  over  the  entire  band.  The  corresponding  aperture  field 
distributions  are  shown  in  Figure  4.  Looking  at  these  aperture  dis¬ 
tributions,'''  one  can  explain  how  the  beam  equalizer  works.  It  goes 
as  follows.  The  aperture  distribution  for  the  E-plane  pattern  Is 
uniform  which  gives  a  familiar  (sin  x)/x  type  of  pattern.  The  H-plane 
pattern  is  due  to  a  cosine-taper  type  of  aperture  distribution  if  there 
is  no  beam  equalizer  present.  Hence,  the  beamwidth  is  larger  than  that 
of  the  E-plane  pattern  for  a  square  waveguide  aperture.  However,  the 
introduction  of  the  beam  equalizer  forces  the  H-plane  aperture  field 
to  vanish  at  the  center  of  the  aperture,  which  makes  the  field  distri¬ 
bution  look  more  uniform.  Hence,  the  main  lobe  of  the  H-plane  pattern 
narrows  to  achieve  the  beam  equalization  effect  in  the  E-  and  H-planes. 
The  pattern  is  mainly  determined  by  the  septum;  the  post  is  present  for 


These  curves  are  generated  using  only  the  first  three  terms  in  the 
modal  series  expansion.  While  the  boundary  conditions  at  the  wave¬ 
guide  walls  are  satisfied  perfectly  by  definition  of  the  modal  func¬ 
tions,  there  appears  to  be  a  small  residue  of  Ex  at  the  septum  in 
the  middle  of  the  waveguide.  This  residue  is  solely  due  to  the 
truncation  of  an  infinite  series  to  a  finite  number  of  terms.  Perfect 
cancellation  can  be  approached  when  more  and  more  terms  are  used  in 
the  series. 
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impedance  matching.  The  computed  reflection  coefficients  over  the 
frequency  band  are  shown  in  Table  2.  If  the  post  is  absent,  the  reflec¬ 
tion  coefficients  will  be  much  larger  than  those  shown  in  Table  2. 


IV.  Conclusions 


A  septum  and  an  Impedance  matching  post  used  as  a  beam  equalizer  in  an 
open-ended  waveguide-feed  for  reflectors  used  in  satellite  communications 
systems  have  been  analyzed  by  using  a  spectral  domain  approach.  The 
computed  radiation  patterns  in  the  E-  and  H-planes,  as  well  as  the 
impedance  match  results,  have  been  presented  in  Che  paper.  The  performance 
of  the  beam  equalizer  over  the  entire  band  of  operating  frequency  has 
been  evaluated.  The  results  indicace  that  the  E-  and  H-plane  principally 
polarized  patterns  are  equalized  extremely  well  over  the  entire  frequency 
band  of  operation  and  chat  the  impedance  matching  is  also  satisfactory. 
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Tab la  1.  BASIS  FUNCTIONS  IN  SPATIAL  AND  SPECTRAL  DOMAINS 


Spatial  Domain 


Spectral  Domain 
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JQ,  J^,  and  J 2  are  Bessel  functions  of  the  first  kind  and  of  order  zero,  one,  and 
two,  correspondingly. 
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TABLE  2 

REFLECTION  COEFFICIENTS  OVER  THE  FREQUENCY  BAND 
d  -  0.1A 


1 

Frequency  (GHz) 

|  Reflection  coefficient  R 

VSWR  -  itiiL 
1-1  R| 

3.75 

0.2729/ -130.25’ 

1.751 

3.35  1 

1 

0.1879/-134.9O0 

1.463 

3.95  ; 

0. 1071/-142. 02“ 

1.240 

4.05 

0.0358/-167.87* 

1.074 

4.15 

0.0442/73.10“ 

1.092 

POST — ^  SEPTUM  —7 


L  -L 


X 

— 

0  z 


LONGITUDINAL-SECTION 


CROSS-SECTION 


a  =1.3  X 
b  =  1.3  X 

c  * 0.25  X  f  =3.95  GHz 

d  IS  VARIABLE  X  *0.076  m 
t  =  0.062  X 


Figure  1.  Geometry  of  the  waveguide  with  septum  and  post. 
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SCATTERING  FROM  FREQUENCY  SELECTIVE  SURFACES 
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ABSTRACT 


In  this  paper,  we  apply  a  novel  technique,  called  the  spectral-iteration 
approach,  for  analyzing  che  problem  of  scattering  from  periodically  perforated 
screens  which  find  useful  applications  as  raiomes,  optical  filters,  artificial 
dielectrics,  and  so  on.  The  formulation  is  carried  out  in  the  spectral  domain 
where  a  sec  of  algebraic  equations  is  obtained  directly  for  the  spectral  coef¬ 
ficients  of  the  aperture  field  distribution  (or  the  induced  currenc  density) 
rather  than  via  an  integral  equation  formulation.  These  equations  are  then 
solved  simultaneously  using  an  Iterative  procedure  developed  in  this  pape. 
that  circumvents  the  need  for  matrix  inversion.  Because  the  matrix  solution 
is  avoided  in  che  spectral  approach,  it  is  capable  of  handling  large  aperture 
sizes  in  a  computationally  efficient  manner.  The  efficiency  of  computation 
results  from  the  use  of  the  FFT  (Fast  Fourier  Transform)  algorithm  which  is 
employed  in  the  derivation  of  the  algebraic  equations  and  in  the  iteration 
procedure.  A  unique  feature  of  the  spectral-iteration  approach  is  chat  it 
has  a  built-in  boundary-condition  check  which  provides  a  reliable  indication 
of  che  accuracy  of  the  solution.  This  paper  also  shows  chat  Che  spectral 

domain  technique  can  be  applied  to  even  a  wider  class  of  geomecries,  e.g., 
the  3tep  discontinuity  in  a  waveguide. 


The  work  was  supported  by  the  Office  of  Naval  Research,  Contract 
N00014-75-C-0293. 


r.  introduction 


Periodic  structures  such  as  arrays  of  conducting  strips  or  periodically 
perforated  screens  which  can  be  either  free-standing  or  printed  on  dielectric 
substrates  (see  Fig.  1)  have  frequency  selective  properties,  and  find  many 
applications  as  artificial  dielectrics,  optical  and  quasi-optical  devices, 
and  dlchrolc  surfaces  for  antenna  reflectors  and  radooes. 

Conventionally,  the  problem  of  electromagnetic  scattering  from  these 
periodic  structures  is  attacked  using  the  mode-matching  procedure  employed 
in  conjunction  with  the  method  of  moments.  A  description  of  this  procedure 
can  be  found  in  a  number  of  papers  on  the  subject  by  Chen  [1],  Lee  [2],  and 
McPhedran  and  Haystre  [3].  Though  this  method  works  quite  well  in  the  low- 
frequency  region,  it  becomes  prohibitively  costly  if  not  impractical  at  the 
high  frequency  region  where  che  aperture  sire  is  one  to  two  wavelengths,  or 
larger,  because  che  matrix  size  required  for  an  accurate  solution  becomes 
prohibitively  large  and  che  numerical  computation  becomes  extremely  time- 
consuming  and  costly.  The  high  frequency  techniques,  e.g.,  GTD,  cannoc  be 
applied  to  circumvent  the  above  difficulty  either,  because  the  complex 
geometrical  configuration  of  the  structure  does  not  lend  Itself  to  the  ray 
formalism  of  GTD.  In  this  paper  we  introduce  a  new  technique  based  on  the 
spectral  domain  approach  which  provides  an  efficient  and  accurate  solution 
to  che  gracing  problems  described  in  this  paper. 

As  a  first  step,  che  new  approach  begins  with  the  formulacion  of  the 
problem  In  terms  of  an  Integral  equation  In  the  transform  domain.  The 
scandard  procedure  for  deriving  the  Integral  equation  for  the  unknown 
aperture  field  (or  the  induced  current)  Is  sccii  followed;  however,  in  the 
transform  domain  the  convolution  form  of  the  integral  eouation  oeccmas 
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FIELD 


SUBSTRATE 


a)  aid*  'rLmt 


?igur«  i.  "r«qu«tiC7  3*i*c  :ivi  Sviriaca. 


an  algebraic  on*.  Furthermore,  because  of  the  periodic  nature  of  the 
structure,  the  transform  naturally  takes  the  form  of  OFT  (discrete  Fourier 
transform)  which  can,  in  cum,  be  efficiently  evaluated  using  the  FT?  (fast 
Fourier  transform)  algorithm.  The  transformed  integral  equacion  is 
subsequently  solved,  using  an  iterative  procedure,  simultaneously  for  the 
aperture  field  and  the  induced  current.  It  is  evident  chat  the  method 
avoids  the  time-consuming  steps  of  evaluating  the  matrix  elements  and 
their  subsequent  inversion.  More  importantly,  the  problem  of  storing  and 
handling  over-sized  matrices  is  circumvented  even  at  high  frequencies, 
where  the  number  of  unknowns  can  exceed  the  figure  2000.  An  added  feature 
of  the  method  is  that  a  built-in  step  in  the  iterative  procedure  provides 
a  convenient  measure  for  the  boundary  condition  check,  a  feature  not 
readily  available  in  conventional  approaches.  Finally,  the  convergence 
of  the  iterative  procedure  is  enhanced  by  combining  it  with  a  variational 
approach  in  a  manner  explained  later.  For  completeness,  we  mention  that 
an  approach  similar  to  the  present  one  has  been  used  by  Mittra  and  So  [ i ] 
in  studying  the  slngle-scatterer  (as  opposed  to  the  periodic  gracing) 
problem.  The  above  authors  have  used  the  GTD  solution  as  the  zeroth  order 
approximate  solution  and  have  also  employed  an  iterative  procedure  to 
generac*  the  final  solution.  However,  to-date  this  procedure  has  not  been 
applied  to  the  gracing  problems  being  considered  in  this  paper. 

In  the  next  section,  we  present  the  formulation  of  the  periodic 
grating  problmm.  In  Section  3,  w*  describe  the  iterative  procedure.  In 
Section  *  w*  illustrate  the  application  of  the  cechnJ qua  to  a  number  of 
practical  geometries.  Finally,  va  demonstrate  in  Section  3  chat  the 
approach  is  usaful  for  a  class  of  closad-reglon  problems,  a.g.,  waveguide 
discontinuities .  A  brief  summary  of  the  paper  is  included  in  Section  6 
and  soma  conclusions  art  prasantsd. 
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II.  FORMULATION 


For  the  salt  a  of  illustrating  the  spectral  approach,  we  consider  tha 
problem  of  a  uniform  plana  wave  scattered  from  a  free-standing  periodically 
perforated  conducting  screen  shown  in  Fig.  2.  However,  the  method  of 
solution  is  easily  and  conveniently  extendable  to  the  case  of  a  screen  on 
a  dielectric  subs crate. 

Due  to  Che  periodicity  of  the  structure,  the  electric  field  on  either 

side  of  the  screen  can  be  expended  in  terms  of  the  Floquet  space  harmonics. 
1 U)  c 

Using  the  eJ  time  convention  (suppressed) ,  we  can  write 
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where  E1  represents  the  incident  field,  (X'r  ,  Y+  )  and  CC”  ,  Y”  )  are  the 

am  am  an  am 

reflection  and  the  transmission  coefficients  of  the  Floquet 's  harmonic 
modes,  respectively  and  p^'s  stand  for  the  Floquec  mode  functions,  given 
by 
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2.  Free-e tending,  periodically  perforated  screen  illuminated 
by  a  plane  wave  propagating  la  it  direction. 
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The  z-componenc  of  che  E-field  can  be  derived  from  (1)  using  che  divergence 
theorem. 

Enforcing  Che  condition  that  che  tangential  field  is  continuous  across 
che  interface,  we  obtain 

X*  -  X~  .  Y*  -  Y~  when  m  j*  0  or  n  A  0 

jqq  am  hi  i  jjjj 

and 

X+  ♦  E1  -  f  ,  Y+  +  E1  -  Y"  . 
oo  x  oo  oo  y  oo 

The  8-field  in  che  region  z  >  0  and  z  <  0  can  be  derived  from 
Equation  (1).  Evaluating  che  3-field  ac  z  ■  0  and  z  -  0+,  subtracting 
the  expression  for  one  from  che  other,  and  asking  use  of  che  requirement 
that  che  tangential  components  be  continuous  across  the  aperture,  one 
arrives  ac  che  equation: 
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where  3*  is  the  incident  3-fleld, 

J  is  che  Induced  current  on  the  surface. 
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la  (2)  we  have  used  the  notation  chat  for  a  function  f(r)  defined  on  the  2*0 
plane,  where  r  is  the  position  vector  on  chat  plane,  the  truncation  operator 
9  is  defined  by 

3(f(r))  -  0  for  r  on  the  conducting  surface 

•  f(r)  for  r  in  the  aperture 

and 

9 (f (r) )  -  f(r)  -  9(f (r) )  . 

The  obvious  identity  9(J)  *  1  and  chat  z  *  [H(z  *  0+)  -  3(s  ■  3  )]  ■  J 
have  also  been  used  in  deriving  (2). 

Unlike  the  integro-dif ferential  equation  in  the  conventional  net hod, 
which  applies  only  in  the  aperture  (or  scrip)  region,  (2)  is  valid  over 
the  entire  surface.  The  price  paid  for  extending  the  equation  to  the  full 
range  is  the  introduction  of  an  extra  unknown  J.  However,  as  we  will  soon 
see,  the  additional  unknown  J  can  be  solved  for  along  with  the  aperture  field 
using  the  iterative  procedure  discussed  in  the  next  section. 
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III.  ITERATIVE  PROCEDURE 


The  summation  Involved  la  (1)  and  (2)  can  be  readily  identified  as  the 
DFT  operation,  let  F  be  the  operator  representing  the  OFT,  and  let  E_ 
represent  the  tangential  electric  field  in  the  transformed  domain. 
Identifying  CX”  ,  l”  )  in  (2)  as  the  Fourier  coefficients  of  E„  ,  and 

3HI  Tjtyi  C 

writing  G  for  the  matrix 


we  can  write  (2)  symbolically  as 

F(C-tc)  -  -H*  +  9(J)  .  (3) 

where  the  subscript  t  Indicates  the  tangential  components,  and  it  is 
understood  that  all  the  quantities  are  evaluated  at  z  ■  0. 

If  the  induced  current  were  available,  the  solution  for  E„  could 
be  immediately  obtained  by  Invoking  (3)  and  by  using 

tc  -  G”1r/1(-5J  +  9(J)) )  .  (1) 

In  practice,  however,  J  is  the  unknown  to  be  solved  for,  together  with  E_ 
and  hence  (i)  cannot  be  used  directly.  Instead  of  using  (i),  a  recursive 
relation  between  the  (n  +■  1) approximate  solution  E^n+^  and  the 
approximation  Is  now  derived  and  the  evo  unknowns  Sc  and  J  ara  solved 

for  simultaneously  using  an  iterative  procedure. 

To  derive  the  recursion  formula  for  E^  ,  we  begin  with  '3),  which 

w 

relates  and  E^a\  and  write 
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(5) 


e(j(a>)  -  p(G*£^a>)  + 


Substituting  (3)  into  (4),  on*  obtain* 

1(0+1)  „  g"L(P‘1(-H*  +  9(P(G-t£a))  +3*))) 


(6) 


?(n) 


Equation  (6)  is  th«  desired  recursive  formula.  3tfor«  inserting  E^  into 
(3),  w*  adjust  its  amplitude  by  multiplying  with  a  scale  factor  K,  computed 
according  to  the  variation  expression 


<S£a),  -H*> 

<S£a),  F(G-f£a))> 


(7) 


where  <f,  g>  -  /  _  f*g*da. 

’  *  aperture  * 

Equation  (7)  is  obtained  by  applying  the  one-term  Gaierlcin'3  method 
to  Equation  (3)  using  as  the  testing  function.  It  is  apparent  chat 

£  «  1  when  f£a^  is  the  exact  solution.  K,  therefore,  also  provides  an 
indication  of  the  accuracy  of  the  nCtl  iterated  result  in  a  weighted- 

average  sense. 

In  the  following- we  procc *d  to  outline  an  iterative  procedure  for 
solving  (6) : 

1. 


Begin  with  an  initial  estimate  .  The  amplitude  of 

w 

-(0) 

E^  Is  to  be  properly  adjusted  -using  the  scale  factor 


2. 


3. 


K  determined  from  Equation  (7) . 

Compute  E^  ,  the  discrete  Fourier  transform  of  E^. 
This  step  can  be  carried  out  efficiently  using  the 
?TT  algorithm. 

-  i(Q) 

Compute  G’£c 
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4.  Obtain  cha  OFT  of  G-E^  using  FFT. 

5.  Subtract  -H^  from  tha  rasult  obtained  from  step  4. 

This  gives  tha  zeroth-order  approximate  solution  J^. 

Generally,  tha  approximate  solution  for  1  obtained  in 

this  step  has  non-zero  values  extending  beyond  the 

conducting  surface.  The  satisfaction  of  the  boundary 

condition  for  the  induced  current  can  thus  be  verified 

cti 

by  checking  how  well  the  n  approximation  for  the 
current  is  confined  to  the  conducting  surface. 

6.  Add  to  obtained  in  the  last  step,  and  take  the 

inverse  DFT  of  the  rasulc  using  FFT. 

--I 

7.  Multiply  G  by  the  result  obtained  from  step  6., 

0.(1) 

obtaining  E^ 

3.  Take  the  inverse  transform  of  E^  to  get  .  The 
exact  solution  for  Eg  should  have  zero  value  on  the 
conducting  surface.  This  criterion  serves  as  a  boundary 
condition  check  for  the  approximate  solution  E^ 
obtained  in  this  step. 

9.  Repeat  the  whole  procedure,  as  necessary,  using  9(E^) 
to  generate  the  next  higher-order  solutions  J*"1'  and 
£c  until  convergence  is  achieved. 

In  the  following  two  sections  several  examples  are  presented  to 


illustrate  the  application  of  the  technique  described  above. 
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IV.  SCATTERING  FROM  GRATINGS  AND  GRIDS 


Lee  us  consider  a  free-standing,  scrip  grating  structure  illuminated 
by  a  normally  Incident  uniform  plane  wave  as  shown  in  Fig.  3.  Let  the 
incident  E-field  be  polarized  parallel  to  the  edges  of  the  scrip  (an 
H-wave) . 

The  formulation  for  this  problem  is  given  by  Equation  (2).  The 
iterative  procedure  discussed  in  Section  III  is  applied  to  solve  for 
the  tangential  aperture  E-field,  E  ,  and  the  induced  current  density, 

J. 

Figure  4a  shows  the  incident  E-fleld  truncated  in  the  aperture,  which 
is  used  as  the  zeroch-order  approximation  for  E  ,  i.e.,  e£°^  -  e(E*). 

The  derived  from  E^  is  shown  in  Fig.  5a.  has  significant 

non-zero  values  extending  into  the  aperture  region.  Thl3  could  be  expected 
because  of  the  crude  initial  estimate  made  for  E„®  .  Figure  5b  shows 
obtained  after  one  iteration.  Observe  the  significant  improvement  achieved 
with  just  a  single  iteration  even  though  the  zeroch-order  approximation  for 
e£^  was  rather  crude.  Higher-order  soiu  cions  for  and  obtained 

via  further  iterations,  are  shown,  respscclvely,  in  Figures  4b  to  4d  and 
Figures  5b  to  5d.  The  rapid  convergence  and  the  accuracy,  vhich  is  verified 
by  che  boundary  condition  check  of  the  solutions,  are  veil-demons  traced 
in  these  figures.  The  induced  current  density  J  also  shows  the  expected 
edge  behavior,  i.e.,  it  becomes  Large  at  the  edges  as  it  should  for  the 
incident  H-wave. 

Figures  6  and  7  show  the  solutions  for  Ec  and  J  obtained  after 
four  iterations  when  che  gratings  ars  illuminated  by  an  obliquely  Incident 
plane  wave  with  an  incident  angle  i  •  30*  and  with  che  H-fleld  polarized 
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Figura  3.  Fr««-scandiag  gracing  tlluolnacad  by  a.  normally  incldenc 


3-wav« . 
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figure  U .  Magnitude  ot  (lie  aperture  E-flelii  distribution  for  the  geonetry 
shown  In  figure  J,  obtained  via  four  consecutive  iterations. 
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gull.-  'J,  obtained  via  four  cunseL-uiive  iterations 


the  induced  current  dentil ty  distribution  on  the 
geowetry  in  Figure  3,  obtained  via  tour 
tcrat ions . 


parallel  co  the  edgaa  of  che  strip  (in  the  E-wave  case).  Again,  the  boundary 

conditions  are  satisfied  extremely  well  by  the  results,  and  the  aoerture 

E-field  also  shows  the  expected  edge  behavior  for  Che  incident  E-wave. 

Next  we  consider  the  scattering  from  a  free-standing  conducting  grid 

illuminated  by  a  normally  incident  plane  wave.  The  geometry  of  the  problem 

2 

is  shown  in  Fig.  3.  The  aperture  area  is  approximately  10X  whereas  the 

2  ^ 

cell  area  is  about  UU\“ .  the  Initial  approximation  for  E_  is  3till  chosen 
co  be  che  truncated  incident  field,  and  che  dominant  component  of  che 
cangencial  aperture  E-flald  is  shown  when  the  incident  E-fleld  is  polarized 
in  che  y-direction. 

For  all  che  computations  in  this  section,  32  terms  in  che  Floquet 
expansion  functions  are  used  in  representing  che  unknown  fields  along 
each  of  che  two  dimensions.  This  leads  co  2^  equivalent  unknowns  to  be 
solved  for.  The  computation  time,  however,  required  for  deriving  che 
solution  is  quite  moderate  (3*6  secs,  of  CPU  time  on  che  CSC  Cyber 
175  System).  Clearly,  any  matrix  method  dealing  with  such  a  large  number 
of  unknowns  will  be  totally  impractical. 


of  t lie  dominant  component 


V.  GENERALIZATION  TO  WO- REGION  PROBLEMS 

la  chis  section,  the  iterative  technique  is  further  generalized  to 
analyze  a  wider  class  of  geometries.  These  geometries  are  characterized 
by  the  feature  chat  they  comprise  a  junction  of  two  or  more  regions  of 
dissimilar  dimensions.  An  open-region  type  example  is  a  corrugated  surface 
which  can  be  thought  of  as  a  junction  of  two  regions,  viz.,  the  infinite 
half-space  and  a  periodic  array  of  short-circuited  waveguides.  For  the 
sake  of  simplicity,  a  closed-region  type  problem  -  a  step  discontinuity 
in  a  parallel-plate  waveguide  -  is  considered  in  this  section.  The  geometry  is 
shown  in  Figure  9.  The  incident  field  is  a  TE  mode  wave.  The  formulation  of 
chis  problem  can  be  found  in  the  literature.  The  integral  equation  is  given  by 
rb  rb  .  , 


fb  -  +  i 

E  (x’)K  (x.x1 )dx'  -  E  (x’)K  (x.x')dx’  -  2H1 

J0  7  J0  7  * 


for  0  <  x  <  b,  z  »  0 


where  Z  is  the  unknown  aperture  E-field, 

y 

is  the  incident  E- field, 

K  (x,x')  ■  -yi—  1  T  sin  —  x  sin  —  x' 

1  u»u  _ ~  m  a  a 

J  a*  ■ 

m 

*  jir  sla  F *  sln  f1  x'  * 

j'j(k2  -  (f^)2)172  if  k2  >  (~)2  , 


tf  k2  <  (fi)2  . 


[  ((f1)2  -  k2)L/2 


! 

*  a  i 


l-j(k2  -  (2L)2)172 

1  3 

if  k2  > 

,3W.  2 

3  ’ 

\ 

!  ( (f2-) 2  -  k2)172 

l.  b 

tf  k2  < 

(f)2 
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Figure  9. 


Step  discontinuity  in  a  parallel-place  waveguide 
with  a  IE  incident  wave. 
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Moce  chac  Equation  (3)  la  defined  in  the  region  0  <  x  <  b.  to  apply  the 
Iterative  technique,  it  has  to  be  extended  to  the  full  range  0  <  x  <  a. 
This  is  achieved  by  introducing  an  extra  unknown  function  J (x) ,  and  the 
extended  equation  takas  the  fora 
a 


» 

E  (x')lf  (x.x)dx’ 
0  7 


fa 

E  (x')K+(x,x')dx’  +  S^H1)  +  i(J(x))  (9) 

0  7  x 


where  for  any  function  f (x) 


for  0  <  x  <  a 


a(f(x))  -  f(x) 

-  o 


if  0  <  x  <  b  , 
if  b  <  x  <  a  , 


and 

9(£(x))  -  £(x)  -  9<£(x)) 


A  recursion  formula  relacing  the  (n  +  l)Ch  order  solution  £^n+'^  to  the 
nCil  solution  can  be  derived  via  a  procedure  similar  to  chat  developed 

in  Section  III.  the  formula  is 


e(oKU- 

y 


E(n)x+  +  H-2ML)  +  e  ( f  E(tl)K  ■ 

0  r  *  J0  7 


E(n):<+ 


9(-2H*))  . 
(10) 


Equaclon  (10)  is  now  solved  using  an  iterative  procedure  similar  to 
that  developed  in  Section  III.  the  integrations  in  (10)  can  again  be 
carried  ouc  using  the  FIT  because  of  the  characteristic  nature  of  the 
kernels.  Figure  10  shows  the  aperture  E-field  distribution  at  the 
discontinuity.  The  Initial  approximation  for  E^  is  taken  to  be  the 
incident  field  truncated  in  the  aperture. 

the  result  is  obtained  in  three  iterations  with  32  expansion  functions 
'used  in  representing  the  unknown  field.  We  note  chac  the  boundary  condition 
on  £^,  viz.,  that  it  vanishes  at  the  edges,  is  satisfied  by  the  iterated 


solution. 
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VI.  CONCLUSION 


la  chis  paper,  an  iterative  technique  applied  in  the  3pectral  domain 
has  been  employed  to  solve  the  problem  of  scattering  from  two-dimensional 
periodic  structures.  Applications  of  the  technique  have  been  illustrated 
using  a  number  of  examples.  This  spectral-iteration  technique  has  been 
found  to  be  very  efficient  for  tvo  reasons:  (i)  it  is  applied  in  the 
transform  domain,  where  only  algebraic  operations  instead  of  convolution 
integrations  are  involved,  (ii)  the  use  of  the  efficient  FFT  algorithm, 
which  is  employed  co  carry  out  the  algebraic  DFT,  is  natural  to  the 
procedure.  The  technique  also  has  a  built-in  boundary  condition  check 
which  verifies  the  accuracy  of  the  iteraced  solution  -  an  important  feature 
not  available  in  conventional  methods.  Furthermore,  the  induced  current  and 
Che  aperture  field  are  solved  for  simultaneously  which  is  another  unique 
advantage  not  available  in  ocher  aechods.  Finally,  as  mentioned  earlier, 
chis  technique  is  particularly  appealing  in  the  high  frequency  analysis 
where  the  moment  method  has  unsurmountabla  difficulties  and  ocher  asymptotic 
techniques  fail.  Before  closing,  we  would  like  to  point  out  that  at  the 
lower  end  of  the  frequency  scale  the  effect  of  mutual  coupling  among 
adjacent  apertures  (or  conducting  strips)  becomes  significant  and 
the  field  distribution  in  the  aperture  begins  to  deviate  substantially 
from  the  incident  field.  Consequently,  it  becomes  difficult  to  choose 
a  good  Initial  approximation  for  the  aperture  field,  particularly  the 
cross-polarized  component,  such  that  a  rapid  convergence  of  the 
Iteration  process  is  ensured.  Fortunately,  however,  an  alternative 
method  also  developed  by  the  authors  [5]  and  called  the  Spectral- 
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Galerkin  approach,  has  been  found  to  be  well-suited  in  the  low 
frequency  region.  This  approach  is  based  on  an  application  of  the 
Galerkin's  method  in  the  spectral  domain.  It  selects  the  proper 
expansion  functions  which  are  analytically  Fourier  transformable 
and  which  satisfy  both  the  boundary  and  edge  conditions  for  the 
unknown  fields.  It  has  been  shown  that  the  required  matrix  size  is 
much  smaller  than  in  the  conventional  schemes  and  the  numerical 
computation  is  very  efficient. 
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APPENDIX  E 


A  SPECTRAL-ITERATION  TECHNIQUE  FOR  ANALYZING  A 
CORRUGATED-SURFACE  TWIST  POLARIZER  FOR  SCANNING  REFLECTOR  ANTENNAS 

R.  Kastner  and  R.  Mittra 

ABSTRACT 

In  this  paper  we  present  an  analysis  of  the  corrugated-surface  twist 
polarizer  which  finds  application  in  the  design  of  scanning  reflector 
antennas.  We  employ  the  spectral-iteration  technique,  a  novel  procedure 
which  combines  the  use  of  the  Fourier  transform  method  with  an  iterative 
procedure.  The  first  step  in  the  spectral-iteration  method  is  the  con¬ 
version  of  the  original  integral  equation  for  the  interface  field  into  a 
form  which  is  suitable  for  iteration  using  a  method  developed  previously 
[ 5 ] — [ 7 ] .  An  important  feature  of  the  technique  is  that  it  takes  advantage 
of  the  DFT  type  of  kernel  of  the  integral  equation  and  evaluates  the 
integral  operators  efficiently  using  the  FFT  algorithm.  Thus,  in  contrast 
to  the  conventional  techniques,  e.g.,  the  moment  method,  the  spectral- 
iteration  approach  requires  no  matrix  inversion  and  is  capable  of  handling 
a  large  number  of  unknowns.  Furthermore,  the  method  has  a  built-in  check 
on  the  satisfaction  of  the  boundary  conditions  at  each  iteration. 


The  work  was  supported  by  the  Office  of  Naval  Research,  Contract 
N00014-75-C-0293. 
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A  SPECTRAL- ITERATION  TECHNIQUE  FOR  ANALYZING  A 
CORRUGATED-SURFACE  TWIST  POLARIZER  FOR  SCANNING  REFLECTOR  ANTENNAS 

I .  INTRODUCTION 

This  paper  presents  an  analysis  of  the  corrugated-surface  twist 
polarizer,  a  device  which  is  used  to  rotate  the  polarization  plane  of 
an  incident  wave  by  90®.  Such  a  polarizer  finds  applications  in  scanned 
reflector  antenna  systems  (Fig.  1)  where  a  rapid  mechanical  scanning 
is  achieved  by  the  movement  of  the  passive,  light-weight  polarizer  [1]. 
The  feed  horn  radiates  a  horizontally  polarized  field,  which  is  reflected 
from  the  "transflector"  radome,  comprising  horizontal  strips,  onto 
the  twist  polarizer.  The  twist  polarizer  is  so  designed  that  upon 
reflection  from  the  polarizer  the  horizontally  polarized  incident  wave 
is  changed  into  a  vertically  polarized  reflected  wave  to  which  the 
radome  is  virtually  transparent.  Thus,  the  radome  also  serves  as  a 
collimating  device  and  the  antenna  is  compact  and  light  weight. 

A  common  type  of  twist  polarizer  comprises  a  set  of  thin  metallic 
strips  or  wires  placed  on  top  of  a  dielectric  substrate,  which  is  approx¬ 
imately  X/4  thick  and  is  backed  by  a  ground  plane  [1],  [2],  [3],  [8]. 

In  this  paper,  we  investigate  the  design  of  an  alternative  configuration 
which  also  can  be  used  for  rotating  the  plane  of  polarization  of  an 
incident  wave,  viz.,  the  corrugated-surface  polarizer.  Its  basic 
mechanism  can  be  explained  very  simply  [4]  by  noting  chat  if  the  period 
of  the  grooves  is  sufficiently  small  compared  to  the  wavelength,  a  plane 
wave  polarized  parallel  to  the  grooves  will  be  reflected  essentially 
from  the  outer  surface  of  the  polarizer,  while  a  vertically  polarized 
wave  will  be  reflected  from  the  bottom  of  the  grooves,  which  results  in 
a  differential  phase  delay  corresponding  to  twice  the  groove  depth. 
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Thus,  if  the  incident  field  were  polarized  at  45°  to  the  direction  of 
the  grooves,  the  horizontal  and  vertical  components  would  have  equal 
amplitudes  and  phases,  and  a  phase  difference  of  180°  will  be  produced 
by  a  groove  depth  of  A/4.  Unfortunately,  the  cell  size  is  typically 
not  small  at  millimeter  waves  and  the  above  result  is  not  sufficiently 
accurate  for  designing  the  polarizer. 

The  spectral-iteration  approach  [5],  [6],  [7],  presented  herein, 
provides  a  method  for  efficient  and  accurate  analysis  of  the  corrugated 
twist  polarizer.  Unlike  conventional  moment-method  and  mode-matching 
techniques,  the  spectral  iteration  technique  requires  no  matrix  inversion 
and,  consequently,  large  cells  can  be  analyzed  without  any  difficulty. 
Furthermore,  the  accuracy  of  the  solution  can  be  conveniently  checked 
at  each  stage  of  the  iteration.  This  is  done  by  determining  how  well 
the  solution  satisfies  the  boundary  conditions  and  by  verifying  the 
conservation  of  energy. 


107 


II.  GENERAL  FORMULATION 


Let  the  corrugated-surface  twist  polarizer,  shown  in  Figure  2,  be 
illuminated  by  a  plane  wave  with  a  wave-number  6  in  the  x-direction. 
Because  of  the  periodic  nature  of  the  geometries,  phase-shift  walls  with 
period  a  can  be  placed  in  the  external  region  (z  >  0)  of  the  structure 
and  the  problem  can  be  reduced  to  that  of  solving  for  the  discontinuity 
between  a  waveguide  with  these  phase-shift  walls  and  the  short-circuited, 
parallel-plate  waveguide  of  width  b  in  the  interior  region  (see  Fig.  3). 
In  two  cases,  viz.,  the  TE  polarization  (parallel  to  grooves)  and  the 
TM  polarization  (perpendicular  to  grooves)  need  to  be  solved  for 
individually  to  obtain  the  differential  phase  shift  between  the  two 
reflection  coefficients.  For  an  incident  wave  polarized  at  45°  to  the 
direction  of  the  grooves,  the  plane  of  polarization  will  be  rotated  by 
90°  if  the  differential  phase  shift  between  the  two  reflection  coeffi¬ 
cients  is  180°. 

The  spectral-iteration  approach  for  this  problem  is  based  on  the 
fact  that  the  Green's  functions  for  both  the  external  and  internal 
regions  are  expressible  in  terms  of  Fourier-type  series.  Thus,  in  each 
of  these  two  individual  regions  the  integral  operator  that  relates 
E-  and  H-fields  is  exactly  invertible.  Denoting  the  operator  for  the 
external  region,  which  is  a  cell  bounded  bv  the  phase-shift  walls,  by 
L_,  and  the  operator  for  the  inner  region  by  L+,  we  have  the  general  form 

L_*  +  L+<1>  -  2Flnc  (1) 
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Inc 

In  which  F  corresponds  to  the  electric  field  incident  from  the  outer 
region  and  \ p  is  the  unknown  magnetic  field  in  the  aperture.  (If,  instead, 
the  unknown  is  the  electric  field  in  the  aperture,  then  F  is  replaced  by 
the  incident  magnetic  field.) 
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III.  TM  POLARIZATION 


In  order  to  illustrate  the  application  of  the  spectral-iteration 
technique,  we  consider  the  case  of  the  twist  polarizer  illuminated  by  a 
TM  polarized  wave.  (The  TE  case  is  discussed  later.)  We  choose  to 
be  the  unknown  electric  field  in  the  aperture  and  derive  the  following 
integral  equation  via  the  usual  formulating  procedure  of  matching  the 
interface  fields. 


j  (  6+tt^— )  x ' 

*(x')e 


-j  (S+n^—)x 

e  a  dx'  +  H(x) 


2HinC(x) 

(2) 


h(x)  is  the  unknown  magnetic  field  outside  the  aperture  where  the 

boundary  condition  on  the  conducting  wall  requires  to  be  identically 

zero.  For  the  sake  of  convenience,  we  have  suppressed  the  factor  -jup 
Inc. 

in  H  and  H  .  Equation  (2)  is  of  the  form  of  (1) . 
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The  operator  corresponding  to  the  external  region 

oo  a 

L.C)  .  I  l 


0  /k2  -  (S  +  n  — 


j(B+A(x’-x) 
e  a  dx'(-) 


(5) 


L)‘ 


is  invertible  in  the  range  0  x  _<  a,  whereas  the  inner  operator 


M->  -!  I 


A 


n“0  .  Cl  ,mir.  2  ..  Jl  ImrrTT 

COt(^k  -  (r— )  d)  /x  -  (r— ) 
ob  ob 


cos  —  x  cos  —  x  dx  (•)(6) 
b  b 


has  the  same  property  in  the  range  0  x  <_  b. 

The  iterative  procedure  for  solving  (2)  will  now  be  given: 

1.  Assume  an  initial  ^(x),  such  that  <|/(x)  *  0  in  the  range 
b  Z.  x  1  a* 

2.  Evaluate  L  (\p)  in  the  range  0  <_  x  <_  a  via  a  two-step  FFT. 

inc 

3.  Compute  2H  -  L  to  obtain  H. 

4.  Evaluate  L+^(H)  using  a  two-step  cosine  FFT  transform 
over  the  portion  0  x  <_  b  to  obtain  the  next  estimate  for 

in  the  aperture. 

5.  Add  zeros  over  b  <  x  <  a  and  return  to  step  2. 


It  is  evident  that  the  procedure  requires  the  evaluation  of  L 


-1 


which  in  turn  requires  multiplication  by  the  factor  cot  Uk}  -  (52-) 2  d). 
This  factor  becomes  small  in  the  neighborhood  of  d  *  X/4  and  for  m  ■  0 
and,  hence,  presents  no  difficulties  in  implementing  the  iteration 
procedure.  Had  we  chosen  to  be  the  electric  field,  a  similar  iteration 
procedure  would  have  required  the  multiplication  by  the  factor 
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tan(*k2  -  (^p)2  d)  and  Che  process  would  have  diverged.  Thus,  it  is  | 

important  to  choose  the  proper  unknown  for  i p  when  using  the  iteration 

procedure. 

The  reflection  coefficient  is  computed  at  each  step  of  the  iteration, 
and  convergence  is  achieved  when  its  absolute  value  tends  to  1.  Since  a 

I 

reflection  coefficient  close  to  +1  is  expected,  a  uniform  distribution  may 
be  chosen  as  an  initial  guess.  It  is  found  that,  even  with  this  initial 
choice,  convergence  is  attained  within  three  or  four  iterations. 

ij 


t  i 


i 
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IV.  TE  POLARIZATION 


In  this  case,  the  magnetic  field  distribution  is  chosen  as  the 
unknown  i (i(x).  We  have  the  following  integral  equation: 


1 

a 


b 

0 


1 


(6  + 


2tt_)2 

a 


il'(x') 


j  (0+n^)x’ 


-  j  (  0+Tt^— ■)  X 

e  a  dx'  +  E(x) 


2E 


inc 


(7) 


with 


2 
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00 


l 

m-1 


j  tan(«42-(^) 2  d) 

— -  -  ^ - (.<«•) 
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sin  - —  x 

D 


.  mr 
sin  - —  x 
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0  <  x  <  b 


E(x)  - 


(8) 


0 


b  <  x  <  a 


and 

E^nC (x)  -  ■ -e-  --BX-  (9) 

/k2  -  82 
o 

inc 

where  a  factor  jwe  has  been  suppressed  from  E  and  E  .  Next  we  redefine 
L_  and  L+  in  a  manner  analogous  to  Equations  (5)  and  (6),  and  outline 
the  iterative  procedure  for  solving  (7)  in  the  following: 

1.  Start  with  an  initial  guess  for  (e.g.,  a  uniform  distribution 
with  a  corresponding  reflection  coefficient  of  -1). 

2.  Compute  L+(t|i)  over  0  <  x  <  b  by  a  two-step  FFT  sine  transform. 

3.  Add  zeros  over  b  <  x  <  a  to  obtain  E  over  the  entire  range 

inc 

0  <_  x  a  and  evaluate  2E  -  E. 

4.  Compute  L  *  of  the  last  result  over  0  <_  x  <_  a  and  obtain  the 
next  estimate  for  <|»(x). 

5.  Go  to  step  2. 


(x) 
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For  the  TE  case  being  considered  here,  step  2  requires  the  evaluation 
of  L+(<iO  ,  as  opposed  to  the  computation  of  L+(^)  in  the  TM  case.  The 
evaluation  of  L_,  in  turn,  requires  a  multiplication  by  tan  -  (^-) ^  d. 
However,  this  does  not  pose  a  problem  since  the  term  m  ■  0,  which  diverges 
in  the  neighborhood  of  d  *  X/4,  is  excluded  in  the  sine  transform.  The 
convergence  in  the  TE  case  is  also  attained  with  a  very  few  iterations. 

The  use  of  the  FFT  procedure,  combined  with  the  small  number  of 
iterations  needed  to  attain  convergence,  result  in  a  considerable  saving 
of  computation  time  as  compared  to  the  conventional  mode-matching  and 
moment-method  techniques  which  require  a  matrix  inversion.  Furthermore, 
the  method  can  handle  unit  cells  which  are  either  large  or  small. 

Naturally,  a  larger  cell  size  necessitates  the  use  of  an  increased  number 
of  samples  so  that  a  typical  sampling  interval  of  about  X / 25  is  maintained, 
and  the  selection  of  a  wider  interval  may  slow  down  the  convergence.  Thus, 
for  a  given  program  that  has  fixed  dimensions  of  FFT  vectors,  i.e.,  a 
fixed  number  of  samples,  there  is  an  effective  limit  to  which  the  frequency 
can  be  increased.  However,  if  the  primary  quantity  of  interest  is  thi 
reflection  coefficient,  the  detailed  information  derived  from  a  fine 
sampling  of  the  aperture  is  not  needed  and  a  smaller  number  of  samples 
is  sufficient.  A  technique  for  improving  the  convergence  for  such  a  case 
is  described  in  the  next  section. 


V.  SPEEDING  UP  THE  CONVERGENCE 


It  should  be  noted  that  often  the  main  difference  between  the  approx¬ 
imate  i(i  derived  after  a  few  iterations  and  the  exact  solution  is  a  complex 
factor.  One  can  partly  compensate  for  this  difference  by  multiplying  tp 
at  every  iteration  with  the  variational  factor  x,  where 


_„inc 

X  -  2F _ > 

<ij/,  Lf > 


(10) 


F  is  defined  in  Equation  (1)  and  L  =  L  +  L+.  The  factor  X,  which  tends 
to  unity  as  p  approaches  the  correct  value,  can  be  regarded  as  another 
measure  of  convergence.  Although  the  use  of  X  is  of  little  help  when 
the  convergence  is  good,  the  incorporation  of  this  variational  factor 
does  improve  the  convergence  quite  significantly  for  an  undersampled  p. 
More  importantly,  it  also  helps  to  achieve  convergence  for  an  otherwise 
divergent  case  when  the  variational  factor  is  not  used.  This  is 
demonstrated  in  Table  1. 


TABLE  1.  CONVERGENCE  RATES,  TE  INCIDENCE,  25  SAMPLES 


Cell  size 
(wavelengths) 

No.  of  iterations 
without  variational  factor 

No.  of  iterations 
with  variational  factor 

.9  and  below 

U 

3 

1.1 

19 

7 

1.3 

Divergent  i 

13 
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VI.  RESULTS 


A  corrugated-surface  twist  polarizer  was  constructed  and  experimentally 
tested.  Its  duty  ratio  b/a  was  0.5,  with  a  *  21  mm,  or  0.6X  at  the  center 
frequency  of  8.5  GHz,  and  d  ■  7.62  mm  or  0.216X.  A  comparison  between  the 
computed  and  measured  axial  ratios  at  various  frequencies  is  presented  in 
Fig.  4  for  an  incidence  angle  of  5°.  Agreement  between  computations  and 
measurements  appears  to  be  very  good,  despite  an  uncertainty  of  about 
tl°  in  the  measurement  setup  and  a  slight  inclination  of  the  plane  of 
incidence  relative  to  the  direction  of  periodicity.  Interestingly, 
when  the  frequency  is  increased  beyond  the  grating-lobe  bound  of  about 
d/X  -  1  (not  shown  in  Fig.  4),  the  computation  predicts  a  sudden  change 
in  the  axial  ratio  and  the  reflection  coefficient,  whose  magnitude  is 
no  longer  equal  to  one.  The  loss  of  reflected  energy  in  the  fundamental 
space  harmonic  is  obviously  accounted  for  by  a  second  space  harmonic, 
or  a  grating  lobe. 
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Figure  4.  Twist  Polarizer  Performance 
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